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RESUMO 


Nesta tese são apresentados e analisados alguns recen- 
tes resultados obtidos na ârea de códigos corretores de erros, re- 
lacionados com decodificação utilizando decisão suave. Tal procedi 
mento & usado com o fim de evitar a degradação no desempenho de 
sistemas codificados que resulta quando na recepção uma quantiza - 
ção abrupta precede a decodificação. 

No capítulo I é apresentado um tratamento introdutório 
sobre códigos lineares, estudando códigos de bloco e códigos convo 
lucionais. Em seguida, no capítulo II, é abordado o conceito de de 
cisão suave e é apresentada uma maneira de utilizá-la, adaptando ! 
técnicas algébricas de decodificação por distância minima a uma me 
dida de distância que permite o uso de informação de verossimilhan 
ça. Alguns procedimentos subôtimos que utilizam decisão suave são 
também apresentados. No capítulo III são estudados dois algoritmos 
ótimos que fazem uso de tais técnicas na decodificação de códigos! 
lineares em canais sem memória.0 desempenho destes algoritmos e 
analisado através de simulação em computador digital, no capítulo" 
IV, para vários: códigos. Os resultados obtidos são comparados com 
procedimentos que empregam decisão abrupta. Finalmente, no capiítu- 
lo V, comentários gerais são feitos discutindo vantagens e restri- 


ções decorrentes do uso destas técnicas, assim como sugestões para 


investigações futuras. 
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ABSTRACT 


In this thesis, some recent results obtained in the area 
of error correcting codes, related to soft decision decoding techni 
ques, are presented and analyzed, Such procedure is used in order to 
avoid the degradation in the performance of encoded systems which 
results when in the reception a hard decision proceedes the decoding. 

In chapter I an introductory treatment of linear codes 
is presented, studying block codes and convolutional codes. Then,in, 
chapter II, the coneept of soft decision is approached, and a way of 
using it is presented, adapting algebraic techniques of minimum dis- 
tance decoding to a distance measure which permits the use of like - 
lihood information. Some suboptimal procedures which use soft deci - 
sion are also presented. 

In chapter III two algorithms with optimal properties 
are studied, which make use of such techniques in the decoding of 
linear codes in memoryless channels. The performance of these algo - 
rithms is analyzed via digital computer simulation, in chapter IV |, 
for several codes. The results obtained are compared with procedures 
which use hard decision. Finally, in chapter V, general comments are 
made discussing the advantages and restrictions implied by the use 
of these techniques, as well as suggestions for further investiga - 


tions. 
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CAPÍTULO I 


INTRODUÇÃO AOS CÓDIGOS LINEARES 


1.1 - INTRODUÇÃO 


Sistemas de comunicação são essencialmente usados para! 
transmissão de informação, e em geral são constituídos por um trans 
missor, um meio de transmissão atravês do qual a informação é envia 
da e um receptor. Ao bloco transmissor estã conectada a fonte de in 
formação a ser transmitida - sinais de voz, sinais de TV, dados de 
saída de computadores, dados de telemetria de naves de exploração ! 
espacial ou dados de telemetria transmitidos de uma instalação ope- 
rada automaticamente para uma estação central de controle. 

A medida que os sinais atravessam o meio de transmissão 
(canal), eles são distorcidos; ruídos e interferências a eles se so 
mam e se torna uma tarefa importante interpretar os sinais quando '! 
finalmente recebidos no extremo receptor. 

O problema básico em um sistema de comunicação se resu- 
me em transmitir a informação atravês do canal ruidoso, de uma ma - 
neira confiável. 

O problema da transmissão confiável de informação tem 
representado um desafio constante para engenheiros e pesquisadores" 


em comunicações. Aqui, a confiabilidade diz respeito à imunidade da 


transmissão a ação do ruído e outras interferências, e não a indeci 


. 


frabilidade de mensagens (criptografia). 

Nas últimas décadas tem havido um grande crescimento na 
utilização de sistemas de comunicação digital, ie, sistemas onde a 
mensagem jã estã na forma digital ou É convertida para este formato, 
como no caso dos sistemas PCM/TDM |33|. Um sistema de comunicação 
digital pode ser descrito como sendo constituído pelos seguintes ' 
elementos: Fonte > Codificador > Canal > Decodificador > Destinatã- 
rio. 

A fonte consiste da fonte original de informação e de um 
codificador fonte, enquanto que o canal consiste de um modulador |, 


um meio de transmissão e um demodulador. As razões para o cresci - 


mento do emprego de sistemas de comunicação digital são bastante co 


nhecidas [11]. Para estes sistemas, uma abordagem eficaz para o pro 
blema da confiabilidade envolve o uso de códigos corretores de er - 
ros. Serã mostrado que a introdução de símbolos adicionais (redun: - 
dantes) ao fluxo de simbolos produzidos pela fonte (dígitos de in - 
formação) permitirã que se execute a deteção e correção de erros às 
custas, evidentemente, de um aumento na complexidade do sistema e ! 
de uma redução na taxa de transmissão de dados (ou de um aumento na 
banda passante exigida). A potencialidade máxima dos códigos corre- 
tores de erros foi estabelecida-em 1948 por C.E.Shannon |35| atra - 
vês do teorema de codificação para um canal ruidoso. Este teorema ' 
estabelece que todo canal tem uma capacidade máxima C e para qual 
quer taxa de transmissão R < C, existem códigos de taxa JR, os 
quais, com decodificação de máxima verossimilhança, tem uma probabi 
lidade de decodificação errônea arbitrariamente pequena. O teorema! 
de Shannon demonstrou a existência de códigos os quais resultam ; 
numa probabilidade de erro na decodificação arbitrariamente peque - 


se do É á ' EN 
na, mas nao indicou como construir tais codigos. 


Devido a presença de sinais ruidosos no canal, como ante 


riormente mencionado, podem ocorrer erros durante a transmissão. Es 
tes erros podem ser esporádicos e independentes, neste caso são cha 
mados erros aleatórios, ou podem ocorrer em "burst" de vários erros 
de cada vez, e então diz-se que o canal tem memória. O último caso" 
não serã objeto de estudo neste trabalho. 

Para combater os possíveis erros que possam ocorrer du- 
rante uma transmissão de dados digitais em um canal de comunicação! 
ruidoso, são empregados côdigos corretores de erros, os quais podem 
ser classificados como códigos lineares ou não lineares. Códigos li 
neares têm os seus dígitos redundantes calculados como soma módulo" 
a |36| dos dígitos de informação, enquanto que códigos não linea - 
res empregam lógica não linear (tais como portas e, ou, não, etc , 
no caso binário). No que segue serão abordados os códigos lineares, 
por sua importância prática e por serem estudados na quase totalida 
de das publicações sobre códigos. 

Entre os códigos lineares existem basicamente duas têc- 
nicas distintas de controle de erros para realização de uma trans me 
missão confiável de dados digitais em canais ruidosos: Os côdigos ' 
de bloco e os códigos de árvore. Dependendo da maneira como a redun 
dância ê adicionada aos dígitos de informação (mensagem), um destes 
dois tipos de código é gerado. Códigos para os quais a redundância" 
em um bloco de dígitos verifica a ocorrência ou não de erros, ape - 
nas naquele bloco particular, são chamados códigos de bloco. Jã os 
códigos onde a redundância em um bloco verifica a existência ou não 
de erros em mais de um bloco são chamados de códigos de árvore. A 
subclasse mais importante destes códigos & a dos códigos convolucio 
nais, os quais são mais simples e de fácil implementação com rela - 
ção a outros tipos de códigos de ârvore. 

Os códigos de bloco e convolucionais são competitivos ' 


em muitas situações. A escolha final de um deles depende de fatores 


como o formato dos dados, o retardo na decodificação, a complexida- 


A 


Ev E 


de do sistema necessário para alcançar uma determinada taxa de er - 
ros, etc. É dada ênfase neste trabalho aos códigos de bloco, muito! 
embora os algoritmos considerados sejam na sua maioria também apli- 
cáveis a códigos convolucionais. 

Nos anos recentes tem havido um grande interesse em es- 
quemas de decodificação utilizando decisão suave, aplicáveis a estes 
dois tipos de códigos, com a finalidade de evitar a degradação no 
desempenho que resulta quando uma quantização abrupta símbolo-por - 
simbolo precede o decodificador. Em sistemas de comunicação via sa- 
têlite, por exemplo, uma diminuição de 1 dB na relação sinal/ruido' 
proprociona uma economia de milhões de dôlares. 

O intuito desta tese é de apresentar, sugerir e anali - 
sar alguns importantes algorítmos de decodificação de códigos linea 
res que façam uso de tais técnicas. Neste capítulo são estudados de 
forma introdutória os códigos lineares, de modo que seja possível '! 
a compreensão das técnicas de decodificação probabilísticas aborda- 


das em capítulos posteriores. 


1.2 - CÓDIGOS DE BLOCO 


1.2.1 - Códigos de Bloco Lineares 


Os códigos de bloco lineares representam sem dúvida a 
parte mais bem desenvolvida da teoria dos códigos corretores de er- 
ros, fato que se deve a sua forte estrutura algébrica, que permite" 
o emprego de ferramentas matemáticas como a teoria das matrizes, e 
a teoria dos campos de Galois, por exemplo. 

o processo de codificação consiste em segmentar uma men 
sagem em blocos de k dígitos e acrescentar a cada bloco n-k di- 


gitos redundantes. Estes n-k dígitos são determinados a partir '! 


dos Xk dígitos de mensagem e destinam-se a deteção e/ou corre» - 


cão-de-erros quespossam vir a ocorrer durante a transmissão. Des- 
4 , 


te modo, o codificador opera de modo independente com cada bloco , 
e a redundância acrescentada serve apenas para detetar e/ou corri- 
gir erros somente no bloco considerado. 

Um fato importante a ser considerado no estudo de códi 
gos de bloco, é que para um comprimento de bloco n e uma taxa '! 
k/n fixadas, o melhor código de bloco linear tem desempenho quase 
tão bom quanto o melhor côdigo de bloco com os mesmos parâmetros ' 
[5 |. A utilização de códigos lineares resulta em simplificação na 
implementação e na análise do desempenho, por isso são os estuda - 
dos quase na totalidade das pesquisas na área. 

Alguns códigos de bloco simples e bastante utilizados, 

tais como códigos de repetição, códigos de peso constante, códigos 
de arranjos, códigos de único digito de paridade e códigos de Ham- 
ming estão descritos nas referências |5 |,|11|. 
Definição 1.2.1 - Um código de bloco linear C(n,k,d) é um conjun 
to de ar ênuplas com elementos em GF(q), chamadas de palavras '! 
código, as quais diferem entre si em pelo menos d posições e for 
mam um subespaço do espaço vetorial Ve de todas as ênuplas defi- 
nido sobre GF(g). 

Deve ser notado conforme definido acima, que o codifi- 
cador se torna proibitivamente complexo para valores grandes de k, 
desde que devem ser armazenadas as a palavras código de n di- 


gitos cada, em uma memória. Isto pode ser evitado lembrando que o 


código é um subespaço vetorial, como serã visto a seguir. 


1.2.2 - Capacidade de Correção 
Dado um côdigo de bloco C(n,k,d), é importante determi 
nar qual a sua capacidade de correção, ie, quantos erros ele & ca- 
.paz de detetar/corrigir, em cada bloco transmitido. Serã mostrado" 


que esta capacidade de correção estã relacionada com o parâmetro d. 


mis 


Portanto, uma importante figura de mérito de um código é a sua dis 
tância mínima d entre duas palavras código, medida em termos de 
distância de Hamming abaixo definida. 


Definição 1.2.2a - O peso de Hamming de uma ênupla u, denotado 
por W(u), é definido cómo sendo o número de componentes não nulas 
de u. 

' 


Assim, por exemplo, se u = (1,0,2,0,0,5,0,1), então 


w(u) = 4. 


Definição 1.2.2b - A distância de Hamming entre duas ênuplas u e 
v denotada por Do (uv), é definida como o número de componentes" 
nas quais estas ênuplas diferem. 

Assim, se u= (1,0,2,0,1,2) e v = EL, bla) a 


dio 


distância de Hamming entre elas é Dy(U,v) 

Se os vetores considerados são binários, da definição! 
de soma módulo 2 é fácil verificar que Dlu,v) = wtuBv), ie 
que a distância entre uev ê exatamente igual ao peso do vetor! 
soma. 

Dado um código linear C(n,k,d), calculando as distân- 
cias entre todos os possíveis pares de palavras código, a menor '! 
distância obtida é chamada de distância mínima do código, denotada 
por d. Portanto, 

d = Min Dy(U,v) Mv Ee € (1=1) 
ufv 


Se ue v são palavras código de um código linear, en 


tão uv tambêm o é, visto que o código é um subespaço vetorial. 
Segue-se deste fato que a distância mínima para um código linear ! 
binário ê igual ao peso minimo dentre as palavras código não nulas, 

Agora, com relação à capacidade de deteção de um códi- 


go de bloco C(n,k,d), torna-se trivial verificar que a ocorrência 


de uma quantidade de erros menor ou igual a d-l resulta em uma ! 


das 229% ênuplas que não são palavras código, o que possibilita 
a deteção de erros. Portanto, um código com distância mínima d é 
capaz de detetar atê (d-1) erros. 

Pode ser facilmente mostrado |36| que na decodificação 
por máximo de verossimilhança com canal BSC, o decodificador esti- 
ma a palavra transmitida como sendo a palavra código mais próxima! 
da palavra recebida no sentido de distância de Hamming. Se o códi- 
go de bloco é usado para correção de erros aleatórios e tem distân 
cia mínima tal que 2t+2 > d > 2t+1, então este decodificador pode 


rã corrigir todas as configurações de erros contendo t ou menos! 


erros que possam ter ocorrido |36 


Assim, um código com distância mínima d tem uma capa 


cidade corretora de t = sl erros. 

Desta forma o problema de encontrar um código com uma! 
dada capacidade de correção, consiste basicamente em garantir uma 
dada distância minima. A teoria da informação proporciona o estabe 
lecimento de cotas superiores e inferiores para d, quando são fi- 
xados um comprimento n e uma taxa k/n, contudo não indica como" 


construir tais códigos |15|. 
1.2.3 - Matriz Geradora de um Código Linear 


Para um subespaço C de V, constituindo um código ! 
linear, é possível encontrar um conjunto de k ênuplas linearmen- 
te independentes, as quais constituem uma base para C. Então, ca- 
da ênupla u e C pode ser expressa como combinação linear dos ve- 
tores (Vj, V,r---+V, ) da base, ie, 

ue c-> um Oy, Og Da Mk (1-2) 
ou seja, 


mo Ve 1 (ve) mp € GF (q) (1-3) 


Um côdigo linear 


V 
n 
dos Xk vetores que 


ser arranjados como 


C(n,k,d) 


constitui um subespaço de 


com dimensão k, de modo que ê possível representá-lo atravês' 


k 
constituem a base. Estes vetores (vph podem 


linhas de uma matriz 


kxn, denominada de 


ma- 
triz geratriz, como mostrado abaixo. 
& nd Mi 
o dE o RO é Ra 
E Saio ad (1-4) 
Ye dy 5 RE 
Se é assumido um bloco de informação m = (my reccmd, 
então a palavra código correspondente é dada por u = mlG]; 
a] k 
u=(mrc cm) = E om V (1=5) 
u 1 Me ts ta UH 
Ye 


ou seja, a palavra código correspondente à mensagem m é uma com- 


binação linear de linhas de [G]. 


Assim, as 


go linear C(n,k,d), 
em uma palavra com n 
& chamada de taxa (ou 


Desde que 


do pela matriz [G], o 


é reduzido enormemente. O codificador pode armazenar somente as 


linhas de [G], ao invês das 


linhas da matriz geradora [G] geram um côdi- 


onde cada bloco de informação é codificado 


dígitos para transmissão. A razão 


R = k/n 


eficiência) do código. 


um código linear é completamente especifica- 


tamanho da memória requerida no codificador" 


k 
q palavras 


k 


código, se ele possuir ! 


elemento lógico capaz de fazer as combinações lineares das linhas! 


armazenadas. 


Exemplo 1.1 - Seja um código de bloco linear (5,3,2) com matriz 


geratriz dada por 


[74 SEDE ni RA 
[Ie] = |v,| = e pe 2/0! RARO 
V, DD sQaa 


A palavra código correspondente à mensagem m = (MOL). 
por exemplo, é encontrada facilmente por 
nbs 


u = (1,0,1) |V,| = 1v, O 07, Iv, EL 0 LÓ 


Es 


É possível codificar cada bloco de informação em pala 
vras código de tal maneira que os Xk primeiros dígitos da pala - 
vra código sejam exatamente os mesmos do bloco de informação, e 
os últimos n-k dígitos sejam os digitos redundantes que são fun 
ção dos dígitos de informação. Um código nesta forma é dito siste 
mático. A redundância é adicionada de maneira a proporcionar capa» 
cidade de proteção à mensagem, e o problema do codificador se re- 
duz a determinar estes dígitos. 

Um código sistemático pode ser descrito por uma ma - 


triz geradora [G] na forma 


[6] =Jo a Ga : (1-6) 


lina 


Onde P43 E GF (q). 

Isto pode ser representado na forma [6] = E Pl, on 
de | é uma matriz identidade kxk e [P] uma matriz kx(n-k) 
com elementos em GF(g). 

Se é considerado um bloco de informação m = (mjreceer 
mM) e uma matriz geradora na forma acima, a palavra código corres 


pondente serã dada por é 
QUO prata (my ve ccem )- [6] & (1-7) 


Efetuando a multiplicação, são obtidas as equações de 


verificação de paridade do código: 


u = 0a 


j (1-8) 
Vas * Pim DO Pr Set, anek 


Os k primeiros dígitos da palavra código são exata - 
mente os dígitos de informação a serem transmitidos, enquanto que 
os últimos n-k dígitos são redundantes (ou de verificação de pa- 
ridade), calculados como funções lineares dos dígitos de informa - 


ção. 


1.2.4 - A Matriz de Verificação de Paridade 


Para cada matriz [G] kxn, existe uma matriz (n-k)xk ! 


[H] de tal forma que o espaço linha de [6] é ortogonal a [H], ie. , 
o produto interno de vetores do espaço linha de [G] com as linhas! 
de [H] é nulo. Esta matriz é denominada de matriz de verificação '! 
de paridade e pode ser escrita sob a forma 
dE 
[H] = [by ++ -Bnp] r onde (Vi) h, & uma ênupla (1-9) 


A 
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Desta forma, se [H] & o espaço nulo de [G], então sê 


verifica a relação 
[6]. [H]' = 0 (1-10) 


Se o código estã na forma sistemática [G] = [IiPJ, en 
tão o espaço linha de [G] &o espaço nulo da matriz [H] dada ' 


T. 
por [H] = [-P kl - 


(1-11) 
Definição 1.2.4 - Se [H] & usada como matriz geratriz de um códi- 
go de bloco, este código é dito ser o dual do código gerado pela 
matriz [G]. 

Seja u = (ujreccuç) um vetor qualquer do espaço li - 
nha de [G]. Então u [H]!= 0 , pois desde que [H] é o espaço nu- 
lo de [G], u.h; = O para 1=1,2, .copnh=k. 

Então o código linear gerado por [G] pode também ser 
descrito de outra maneira: 

"u é palavra código se e somente se utm)"= 0, onde [H] 
é a matriz de verificação de paridade do código linear gerado por 
[G]". 

Exemplo 1.2 - A matriz [H] para o código de bloco (5,3,2) introdu 
zido no exemplo 1.1 pode ser determinada pela equação (1-11) resul- 


tando em 


E 
" 


A énupla u = (1,0,1,1,0) resulta em ulH]'= 0. 
Conclui-se, portanto ,que um código linear é unicamente" 
especificado pela matriz geradora ou pela matriz de verificação de 


paridade. 
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Na construção de um código linear, a matriz [P] deve! 
ser escolhida de modo que o código tenha as propriedades de corre - 
ção de erros desejadas. 

Para um código com matriz geradora [G] = [Le Pl, as 
equações de verificação de paridade podem também ser obtidas direta 


mente da matriz [H] = ei RE dis como mostrado abaixo. Seja 


[e 
1 


(uy rece ru) uma palavra código correspondente a mensagem 


do RR UR a 


I3 
H 


(my... Mm), o que implica em u = M, 


Desde que u[H)” = 0, segue-se que 
U43 = P45 u, 6) RIR O Ps ui «de, 


Ea Pim 6) +++ (O Pes para Jj=1,2,...,n-k (1=12) 


que são as equações já encontradas em (1-8). 
1.2.5 - Sindrome 


Suponha que uma palavra código u de um código de blo- 
co linear C(n,k,d) com matriz geradora [G] e matriz de verifica 
ção de paridade [H] & transmitida através de um canal ruidoso. No 
receptor, uma ênupla r é recebida, a qual pode diferir de u de- 
vido ao ruído adicionado durante a transmissão. A tarefa do decodi- 


ficador é recuperar my) a partir de p. 


Definição 1.2.5 - O vetor de n-k componentes s ser” é chamado 


de sindrome da ênupla r. 


Como visto na seção anterior, yu é uma palavra código 


se e somente se sua sindrome é nula. 


Portanto, a sindrome de um vetor recebido na saída do 
canal pode ser usada para deteção/correção de erros. 


O processo de decodificação envolve uma decisão sobre 


qual foi a palavra código transmitida. O arranjo padrão |24] suge- 


A 


ug 


: : n 
re uma maneira de como fazer isto. No espaço Vp! àS 2 ênuplas ! 


e - k á nd à 
sao distribuidas em 2 conjuntos distintos, de modo que cada um ! 


deles contenha apenas uma palavra côdigo. Todos os elementos de um 
dado conjunto possuem a mesma sindrome, e dois conjuntos distintos 
são associados a síndromes diferentes. Desta forma é estabelecida '! 
uma correspondência um a um entre uma configuração de erros (coset 
leader) e uma sindrome. Neste caso, o decodificador corrige erros ! 
quando a configuração de erros do canal é um "coset leader". Estes 


fatos implicam em um procedimento geral para decodificação de códi- 


gos de bloco lineares, denominado de busca sistemática |30]. 


1.2.6 - Códigos Cíclicos Binários 


Uma considerável parte das pesquisas em côdigos de blo- 
co estão concentradas numa subclasse conhecida como códigos cícli - 
cos, introduzidos por Prange |26 | em 1957. 

Estes códigos são também os mais importantes do ponto 
de vista de aplicações em engenharia. Tudo isto & devido a sua for- 
te estrutura matemática que permite uma considerável simplificação! 
na implementação dos sistemas. Os procedimentos de decodificação de 
códigos de blocos lineares são tambêm aplicáveis a côdigos cíclicos, 
todavia, propriedades algébricas decorrentes da estrutura cíclica ! 
permitem simplificações importantes, visto que a codificação e (o) 
cálculo da sindrome podem ser facilmente realizados empregando re - 
gistros a deslocamento com conecções de realimentação. Sua estrutu- 
ra inerentemente algébrica também torna possível encontrar vários 
mêtodos de decodificação simples e eficientes. Entre os procedimen- 
tos de decodificação mais importantes para este tipo de código, são 
conhecidos DECODIFICAÇÃO DE MEGGITT |36 |, DECODIFICAÇÃO POR LÓGICA! 
MAJORITÁRIA |28| e o método conhecido como "ERROR-TRAPPING" 24] 


os quais são de implementação simples e bastante atrativos. 


ig 


Definição 1.2.6 - Um código de bloco é cíclico quando uma permuta - 
ção cíclica, aplicada a qualquer de suas palavras, resulta ainda em 


uma palavra código, ie, se v = (vor env ) é uma palavra cô- 


n-1 


: e i do ; - - ' 
i ent = : SE Medrado é tambem &, consi- 
digo, entao v RD VorV4r . RE eo o e, co 
derando os indices reduzidos módulo mn. 
É possível tratar as componentes do vetor vw como coefi 
cientes de um polinômio de grau sempre menor ou igual a n-1, atra- 


vês de uma correspondência um a um: 


n: 


W = (vgreccrV no) <=>v (x) = VgtVjXte cet 4X (1=13) 


Fazendo uso de propriedades de campos finitos, prova-se' 
I36 |gue todas as palavras de um código cíclico (n,k,d) são múlti - 
plas de um polinômio g(X), de grau nx, bem definido; e reciproca 
mente que todo polinômio de grau igual ou menor que n-1, e que se- 
ja divisível por g(X), & uma palavra código. 

O polinômio g(X) é chamado polinômio gerador do código 


cíclico e é fator de xº41, ie, 
3 n 
g0)h() = Cr, a TsQnl) = 0 mod XAL. 


A representação de um código cíclico pode ser feita atra 
vês do seu polinômio gerador g(X), da mesma forma que a matriz [G] 
representa um código de bloco linear. 

Conforme mencionado, as palavras código são múltiplas do 


polinômio g(X). Desta forma, os polinômios g(X), EE é PR co 


g(x) 
são palavras código e como são linearmente independentes, podem ! 
ser usadas como base para formar a matriz geradora do código cicli- 


co, como mostrado a seguir. 
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x" Lg (x) 


[]=| x'g(x) 


g(x) 


Para fins de codificação, a propriedade do deslocamento 
cíclico permite uma implementação sequencial da matriz [G]. 

Representando a mensagem a ser codificada por um polinô 
1 


mio m(X) de grau menor ou igual a k-l, um codificador na forma 


sistemática gera a palavra código de acordo com 


v(x) = ct) + xPÉ m(x) (1-14) 


onde c(x) & o resto da divisão de nei m(x) por g(X), 


c(x) = rem(xÉ m(x)/g(x)). (1-15) 


A implementação. do codificador pode ser feita através de 
registro a deslocamento com n-k estágios, pois multiplicações e 
divisões de polinômios com coeficientes em GF(2) podem ser efetua 
das facilmente com o auxílio de registradores a deslocamento |11]. 

Um codificador com k estágios ao invês de n-k está - 
gios pode também ser implementado |30|, levando em consideração o 
polinômio de verificação de paridade h(X). 

A sindrome pode ser facilmente determinada, conforme se 
rã mostrado a seguir. De um modo geral, a palavra recebida pode ser 
representada pela expressão r(X) = v(X) + e(X), onde e(x) é um 
polinômio correspondente à configuração de erros introduzida no ca- 
nal. A síndrome pode ser determinada como o resto da divisão entre 


o polinômio recebido e o polinômio gerador, ie, 


vix) ; ex), intao 


s(x) = remír(x)/g(X)) = rem(—— 
g(x) g(x) 


. 
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Como V(X) = m(X) 9(X), segue-se que 
s(x) = rem (e(x)/g(x)) (1-17) 


Caso S(xX) seja zero, o decodificador aceita a palavra 
recebida como pertencente ao código. No caso contrário, ie, S(X)£Z0, 
considera que ocorreram erros durante a transmissão. Desta maneira" 
é fácil perceber a simplicidade de um circuito decodificador para 


deteção de erros quando são empregados códigos cíclicos. 


1.2.7 - Códigos de Hamming e Códigos B.C.H 


Os códigos de Hamming, introduzidos em 1950 ELG, foram 
os primeiros códigos não trivais propostos para corrigir erros. São 
códigos lineares que corrigem um erro por bloco e têm distância mi- 
nima igual a 3. Para cada inteiro c>1, existe um código de Hamming 


com os seguintes parâmetros: 


A condição n = 2.1 & escolhida de modo a assegurar 
disponibilidade de redundância suficiente para verificar a ocorrên- 
cia de um erro por palavra, pois o número de sindromes não nulas '! 
2º.1 resulta exatamente igual ao número de posições n em que (o) 


erro pode estar localizado. Isto significa que os códigos de Hamming 


são códigos perfeitos conforme a definição abaixo. 


Definição 1.2.7 - Um código de bloco (n,k,d) corretor de t er - 
ros e definido em GF(q) é perfeito se e somente se 
n-k 


does 
Rea Ea 


= po 


Com excessão dos códigos de Hamming, o código de Golay' 
(23,12,7) com t=3 e o código ternário (116.7) com t=3, não exis 
tem outros códigos perfeitos não triviais. 

É interessante observar que os códigos de Hamming são 
também códigos cíclicos e que a taxa (eficiência) é dada por : 
R=1- a o consequentemente são códigos cuja eficiência tende" 
al assintGkicansino. Os códigos (7,4,3), (15,11,3) e (31,26,3) se- 
rão utilizados nas análises realizadas nos capítulos seguintes. 

Com relação aos códigos BCH, estes códigos foram cria - 
dos independentemente por Hocquenghen (1959) e por Bose-Chandhuri ' 
(1960); são códigos cíclicos, e representam a classe mais importan- 
te de códigos de bloco com algoritmosalgébricos de decodificação . 


Para quaisquer inteiros positivos m e t (t<22-1) existe um códi 


go BCH com os seguintes parâmetros: 


O teorema fundamental dos códigos BCH é a seguir enun - 


ciado. 


Teorema 1.1 - O código BCH cujo polinômio gerador tem d-1 raizes! 
consecutivas q, É, É, ia tem distância mínima pelo menos d. 


pocos 


Para um tratamento completo destes códigos o leitor de- 


ve consultar as referências | 3| e |26 


1.3 - CÓDIGOS CONVOLUCIONAIS 
1.3.1 - Considerações Gerais 


Na classe de códigos denominada de códigos convolucio — 
nais, a redundância adicionada a cada bloco é função de mais de um 


bloco de informação, de modo que estes não são tratados de forma in 


. 


Sig 


dependente como nos códigos de bloco. Estes códigos foram primeira- 
mente introduzidos |26 | por P.Elias (1955),e podem ser classifica - 
dos como fixos ou variantes no tempo. Entretanto, a grande maioria" 
dos códigos convolucionais empregados é invariante no tempo, de for 
ma que apenas estes serão abordados aqui. 

Na decodificação de códigos convolucionais, vários sub- 
blocos sucessivos recebidos são processados. Idealmente, toda se -— 
quência transmitida deveria ser empregada na decodificação, o que 
não é feito devido a limitações práticas e por introduzir um retar- 
do muito grande na decodificação. Todavia, a degradação resultante" 
desta simplificação ê inteiramente aceitável. 

Num sentido restrito, códigos convolucionais podem ser 
vistos como uma classe especial de códigos de bloco lineares, por - 
rêm, uma observação mais cuidadosa revela que a estrutura convolu - 
cional proporciona a um código linear propriedades superiores que 
melhoram o desempenho. 

A estrutura de um código convolucional pode ser exibida 
por meio de diversas representações, visto que um codificador convo 
lucional fixo pode ser considerado como uma máquina sequencial de 


estados finitos, invariante no tempo. 
1.3.2 - O Codificador Convolucional 


Na entrada do codificador é alimentado um bloco com k 
simbolos de informação (elementos de GF(q)), e um bloco com n>k 
ptnbidds código é gerado e aparece na saída do referido codificador. 
Os simbolos de saída são combinações lineares em GF(q) dos símbolos 
de entrada dos N-1 blocos precedentes. Deste modo, os n dígitos 
de saída dependem não somente dos k dígitos de mensagem de um mes 


mo sub-bloco , mas também dos N-1 sub-blocos de mensagem anterio- 


res, 


E 


Um código assim gerado é dito ser um código convolucio- 
nal (n,k,N), com comprimento de restrição no codificador de JN t 
blocos. A taxa assintótica para este código é R = k/n, e em geral, 
k e n são inteiros pequenos. 

A formulação matricial para estes códigos é descrita a 
seguir: 

Inicialmente é considerado um conjunto de k.(n-k) ve- 


tores com N digitos binários cada, denominados SUBGERADORES : 
(1,9) =(gg (1,5) r gy (115) recerInog 1d) (1-18) 


Para is DD ssa e j = DS on 


A partir deste conjunto, forma-se uma matriz com k 1i 


nhas, 
To (1) 
ê É zeros  .. o 
: = [1 Pq0P,0P,...0 P4 40 > 00] (U=19) 
Io (k) 
Onde 1. tem dimensão Kkxk 
0 tem dimensão Kkxk 
P, tem dimensão kx(n-k), O<t<n-1 
e 9 pt1,1) 9 (1,2)... .0 (Lyn-k) 
[Pp] = : : : (1-20) 


9 1) gp(k,2)...9 plkynck) 


Assim, cada linha g (i) é um vetor semi-infinito com 
componentes não nulas confinadas aos primeiros nN dígitos. 

O vetor g(i) contendo somente as primeiras nN posi- 
ções do vetor gli) é dito um GERADOR. 

Um código convolucional fixo linear pode ser representa 


do por uma matriz geradora sob a forma: 


De 


zeros, 
LI Pç 0 P, Eis et a > 
zeros 
G. = Lo E) O ond mi nie dam 0 Py-1 > (1=21) 
zeros 
Ie Por CRER e DV 0 Pu as 
y 
co 


Alguns exemplos apresentados esclarecem como determinar 


a matriz G, a partir dos subgeradores. 


Exemplo 1.3 - Considere um código convolucional sistemático (2,1,4) 
o qual é gerado de acordo com o subgerador g(1,1) = (1101). 
O gerador é g(l) = (11010001) ea matriz gera- 


dora ê determinada como sendo 


8+e 


Portanto, para uma mensagem a ser codificada 


= Cl 6 ra a SER À sequência codificada é encontrada por 


H 
3 
q 


je, CAL GOO O O ada 


A Ea (o " 


Exemplo 1.4 - Considere agora o código convolucional (4,3,3) com os 
seguintes subgeradores: 


gttiy = (0 1) 
gt2,d) = (1002) 


g(3,1) (1/20) 


Os geradores são determinados como sendo: 


o 


(1 Q0 a q0. 0.2 Do q) 


g(1) 
g(2) 
gt) = (00 11 0:00 1 0000) 


que 1 Gr 0000 070.0. 1) 


e a matriz geradora serã dada por 


TO Qro pranosi” “o Dios 
SEO 002080==0 0=0:1 zeros 


00 1:1 000:1 60 0:0 


Sos É 100:1 000:1 00 0:1 
010:1 00 o:o 00 Us 
00 lil 000:1 00 0:0 
: 
Deve ser observado que PQ bu 0 9 RO A DO E pj = (oo e 


P; = (110). 

Um côdigo convolucional [n,k,N) pode também ser especi 
ficado pela sua matriz de verificação de paridade H,, tal como um 
côdigo de bloco. Novamente ê verificada a relação G, Ho = 0, donde 
segue-se que c é uma sequência código se e somente se c H =0. 

Para um código convolucional sistemático, ou seja, com 


matriz geradora da forma descrita em (1-21), a matriz de verifica - 


ção de paridade & dada por 


E 4 
Pq Lck 
T ei 
Po O o Sã 
a Au mM 
Pp 0% Mas 
e ' (1-22) 
pg PE 
at 
tao Pe pese PM RR DE ER 
| 
Sp TE JE 
Pyiy- Praia did DO UR 
+ -> o 


Exemplo 1.5 - 


para o código (4,3,3) descrito no exemplo 1.4 é determinada abaixo. 


Como um exemplo, a matriz de verificação de paridade 


Medo Ao, E 
1.010 a 
H e O .0L a LL AI 
1 1 00 "0 10 pe TR E 
11690 :90 106 Ya 
+ -> co 
1.3.3 - Outra Representação da Matriz Geradora 

É possivel tambêm representar a matriz G, numa | outra 


forma que facilita o projeto dos circuitos codificadores, como serã 


mostrado adiante. 


Para um código convolucional. (n,k,N), a matriz G pos 
de ser escrita sob a forma: 
Gg E aRdo es dc 
Ia Me “* Iy1 
Gu (1-23) 
Ig GB errrerees 9n-1 
+ -> o 

| o 

Onde  9;r j = 0,1,... N-1 são matrizes kxn. 


Se o código ê sistemático, comparando (1-21) e (1-23) , 


é possível observar que existem relações da forma 


dq * Tp Po + CEA O ey 4 Pa - 


A matriz de verificação de paridade também pode ser es- 


. 


=23- 


crita na forma abaixo, 


hoy 
h by 
go 4 ba he & (1-24) 


h h araras dA h 


> o 
+ 
oo 


E ST “st “st 
Onde hy=Pg 1h.pr hy5Pj O, ho=P> 0h 000 BgogPyoy 


Os k geradores podem ser facilmente obtidos a partir! 


da relação 


941) 


Ê FO qua rede 
g(k) 


(1-25) 


Os exemplos apresentados abaixo mostram como esta repre 


sentação é refletida na implementação do circuito codificador. 


Exemplo 1.6 - Um circuito codificador para o código convolucional! 
(2,1,3) especificado através de 99" e 9, abaixo,ê mostrado na 


figura Lol. 


9 & (2,0) => gerador glyr= (11 10 11) 
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Figura 1.1 - CODIFICADOR PARA O CÓDIGO CONVOLUCIONAL 
(2173). 


Exemplo 1.7 - Um circuito codificador para o código convolucional'! 
não sistemático (3,2,2) especificado atravês das matrizes 9% € 9' 


abaixo, estã apresentado na figura 1.2. 


e 
ga = 
A RE g(1) =(111101) 
=> geradores 
ba g(2) =t0 10110) 
9 & 


Figura 1.2 - CODIFICADOR PARA O CÓDIGO CONVOLUCIONAL 
Rg E7A Ti 


Portanto, para implementação do codificador de um códi- 
go convolucional (n,k,N) são necessários k.(N-1) estágios de re 


gistro a deslocamento com ligações determinadas pelos geradores. 


. 


espe 


1.3.4 - Diagrama em Árvore, Treliça e Diagrama de Estado 


Uma das maneiras de representar a sequência de saida re 
sultante da codificação de uma dada sequência de entrada para um cô 
digo convolucional fixo, é por meio de uma representação da árvore" 
associada ao codificador. 

Os bits de entrada do codificador são indicados pelo ca 
minho seguido na árvore, enquanto que os bits de saída são indica - 
dos nos ramos. Um zero de entrada é representado pelo ramo superior 
de uma bifurcação, enquanto que o um é representado pelo ramo infe - 
tior. 

O diagrama em árvore para um código convolucional não 


sistemático (4,1,3) é apresentado na figura 1.3 abaixo. 


0000 


Figura.1.3 - DIAGRAMA EM ÁRVORE PARA O CÓDIGO 


CONVOLUCIONAL (4,1,3). 


A 
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trada do codificador produz uma sequência de saída , 

Observando-se o diagrama de árvore, nota-se que sua es- 
trutura torna-se repetitiva apôs um certo número de ramos, o que 
estã relacionado com o fato do codificador convolucional ser uma mã 
quina com estados finitos. Isto torna possível fazer uma modifica - 
ção que consiste em juntar os nôs indicados por uma mesma letra, fa 
zendo com que o diagrama em árvore evolua para um diagrama de treli 
ça. 

A treliça para o código (4,1,3) representado na figura 


1.3 é facilmente obtida, resultando no diagrama mostrado abaixo. 


1001 1001 1001 1001 


Figura 1.4 — DIAGRAMA DE TRELIÇA PARA CÓDIGO 
CONVOLUCIONAL (4,1,3). 


Pode ainda ser observado que ê possível fazer uma repre 
sentação em diagrama de estado de um codificador convolucional. A 
partir da treliça é possível construir este diagrama, que represen- 


ta perfeitamente o codificador. O diagrama para o exemplo em ques - 
tão é mostrado na figura 1.5. a 


mada 


c=[10] 


b=[01] 


1001 


Figura 1.5 - DIAGRAMA DE ESTADO PARA O CÓDIGO 
CONVOLUCIONAL (4,1,3). 


Utilizando este diagrama é possível a determinação do 


conjunto de pesos dos caminhos do código de uma forma relativamente 


simples |37|. 
A determinação da distância livre para um código convo- 
lucional, a qual é obtida a partir do conjunto de distâncias entre" 
o caminho todo zero e os demais caminhos que dele divergem, também" 


pode ser feita atravês do diagrama de estado |30|- 


1.4 - DECODIFICAÇÃO DE CÓDIGOS CONVOLUCIONAES 


Vários procedimentos são usuais na decodificação de có- 
digos convolucionais. Entre as técnicas mais comumente utilizadas ' 
podem ser citadas a DECODIFICAÇÃO SEQUENCIAL |40|, (algoritmo de 
FANO/WOZENCRAFT) e o ALGORÍTMO DE VITERBI |37|, entre outras. 

Para alguns códigos, a decodificação também pode ser 


realizada utilizando LÓGICA MAJORITÁRIA |36|. 


O leitor interessado nestes procedimentos deve se diri- 
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gir à literatura apropriada |1 |, |28|,|40|. 

Ainda com relação a decodificação de códigos convolucio 
nais, serão apresentados no capítulo III, duas maneiras de utilizar 
técnicas de decisão suave. Inicialmente, um decodificador serã apre 
sentado, o qual utiliza decisão suave e ê ótimo no sentido de que 
minimiza a probabilidade de erro por bit transmitido. É também co - 
nhecido que o algoritmo de Viterbi para decodificação de códigos 
convolucionais maximiza a probabilidade a posteriori das sequências 


código em canais sem memória. 


CAPÍTULO II 


DECODIFICAÇÃO USANDO DECISÃO SUAVE 


2.1 - INTRODUÇÃO 


Como jã discutido no capítulo anterior, sabe-se que na 
prática,empregam-se largamente algoritmos de decodificação basea - 
dos em propriedades estruturais de códigos algébricos, o que permi- 
te frequentemente uma implementação relativamente simples em termos 
de circuitos digitais. Deve ser observado, entretanto, que estes 
procedimentos apresentam uma séria deficiência, uma vez que assumem 
um modelo de canal no qual em cada período de tempo um dos q pos- 
síveis sinais do alfabeto ê transmitido e um deles é recebido. 

Este modelo implica em um receptor que efetua uma deci- 
são abrupta, baseada no sinal recebido, escolhendo das q letras 
do alfabeto do transmissor aquela que foi a mais provavelmente ' 
transmitida. Tal receptor descarta,portanto, toda informação acerca" 
da confiabilidade desta escolha, bem como acerca das probabilidades 
das outras letras que não a escolhida. Este fato torna bastante cla 
ro que o desempenho dos sistemas codificados pode ser melhorado pe- 
lo uso da informação probabilística associada à palavra recebida. 

Numa recepção que utiliza decisão abrupta, esta informa 
ção é perdida na quantização que precede a decodificação. Um recep- 


tor que emprega decisão suave é aquele que tenta reter toda (ou par 


te) da informação probabilística para uso apropriado pelo decodifi- 


cador. ; A 
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Em geral, os códigos corretores de erros empregados em 
receptores práticos são côdigos lineares binários. A grande parte '! 
dos procedimentos existentes para decodificação deste tipo de códi- 
go assumem a existência de canais com saída binária, ou seja, empre 
gam um quantizador antes da decodificação para extrair a informação 
na forma digital. A maioria dos detetores existentes constam basica 
mente de um demodulador analôgico, seguido por um decodificador di- 
gital que opera sobre os dígitos binários produzidos pelo demodula- 
dor. O bloco demodulador atua como um quantizador que apresenta um 
limiar de decisão, resultando em duas saídas distintas. Isto signi- 
fica que de acordo com o valor da amostra que é colhida (acima ou 
abaixo do limiar de referência) ,um dígito binário correspondente (1 
ou 0) é entregue ao decodificador. É desta forma que a informação ' 
probabilística, potencialmente útil, entregue ao demodulador & in - 
teiramente destruida, fato que reduz substancialmente a eficiência" 
do sistema de comunicação. Um processo de deteção (demodulação e 


decodificação) com decisão suave utiliza a informação de natureza ' 


probabilística associada. à mensagem recebida, de modo a obter uma ' 
melhor estimativa da mesma. A idêia é fornecer ao decodificador in- 
formação sobre a confiabilidade dos dígitos recebidos, evitando as- 
sim a degradação que resulta quando se utiliza antes da decodifica- 
ção uma quantização de apenas duas regiões. 

A maioria da perda de informação pode ser recuperada se 
a saída do demodulador & quantizada em 2º regiões, 2º-1 em cada 
lado do limiar de referência.Então,para cada dígito de entrada modu 
lado,o circuito demodulador fornece um dígito de decisão (indicando 
se a saída estã acima ou abaixo do limiar),e Q-1 dígitos de confia- 
bilidade (indicando quão distante a saída estã do limiar de decisão) 
Os Q dígitos de decisão suave constituem então a saída total,ao in- 


vês de somente um digito,como no caso de decisão abrupta. 


Na figura 2.1 é mostrado um exemplo para Q=2. 


. 


dies 


1 0 1 O —- Diígito de decisão , 
É 0 0 1 — Dígito de confiança|saida 


Limiar de 
confiabilidade Limiar de decisão 


Instantes de 
amostragem 


Figura 2.1 - DECISÃO SUAVE - QUANTIZAÇÃO EM 4 
REGIÕES DE CONFIABILIDADE. 


A mais simples forma de decisão suave & o apagamento,ou 
deteção por zona nula. Para um estudo mais detalhado, as referênci- 
as |4| e |5 | são apropriadas. Neste método, a região próxima (ime 
diatamente acima e imediatamente abaixo) do limiar de decisão é 


chamada de zona nula. Saídas caindo dentro desta região são entre 


gues ao decodificador rotuladas como apagamento E, como observado 


no exemplo da figura 2.2. 


Limiar de — Zona nula 


decisao 


Figura 2.2 - DETEÇÃO POR ZONA NULA. 
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Quando o decodificador utiliza decisão abrupta, é assu 
mido um modelo de canal BSC, enquanto que o uso de apagamento re - 
sulta em um canal BEC. Os algoritmos algébricos de decodificação ! 
podem ser modificados para manipularem eficientemente os apagamen- 
tos. 

O decodificador utilizando deteção por zona nula tem 
agora algum conhecimento de onde os erros provavelmente ocorreram, 
e pode decodificar de acordo com este fato, ou seja, os d-l1 er - 
ros detetados são localizados atravês da informação de confiabili- 
dade. 

Desta maneira, a potência de correção de erros de um 
código pode ser aproximadamente dobrada, desde que um código com 
distância mínima d pode corrigir d-1 apagamentos, mas somente ' 


d= 
| 3 | erros. 


A utilização de mais de duas regiões de quantização é 


uma tentativa de diminuir a perda de informação que resulta quando 
uma decisão abrupta ê realizada. 

Portanto, o detetor ótimo, ie, aquele que no processo! 
de decisão retêm toda a informação contida no sinal recebido, se - 
ria aquele que entregaria ao decodificador o valor exato das amos- 
tras colhidas. Mais tarde, no capítulo III, estes detetores serão" 
abordados. 

Uma maneira de se utilizar os dígitos de confiabilida- 
de no processo de decodificação é atravês do conceito de distância 
suave. 

A distância suave d entre duas palavras código pode 
ser calculada de acordo com o procedimento descrito abaixo: 

1) - invertendo os dígitos de confiança da palavra se seu dígito '! 
de informação correspondente ê zero. 


2) - Somando módulo 2 as duas palavras modificadas de acordo com o 


Creatas 


passo 1. 

3) - Dividindo a sequência resultante em n subsegquências de q ! 
digitos cada. 

4) - Interpretando cada uma das subseguências de Q dígitos como um 
número binário e convertendo para decimal. 

5) - Somando os números obtidos em 4, resultando em um decimal to - 
tal. 

Assim, d, entre 000 111 e 011 011 é obti 


da por 


O uso de distância suave no processo de decodificação ! 
serã ilustrado a seguir. 

No exemplo abaixo & assumido Q = 3 e um código com 
exatamente duas palavras, 0 ou 1, para o qual d=l1 e R-1. 

A figura 2.3 mostra as possíveis saídas obtidas apôs o 


demodulador . 


limiar de confiabilidade 


limiar de decisão 


digitos dígitos 
de deci de confiança 
são À 
Figura 2.3 - SAÍDAS DO DEMODULADOR EMPREGANDO 
DECISÃO SUAVE. Ê 


ep 


É admitido que a palavra 1 & transmitida, e que devido" 
ao ruído do canal, 0 00 & a saída do demodulador. As duas possibi 
lidades de plena confiança para a saída do demodulador são 011e 
111, de maneira que a distância de decisão suave d+ para cada 


caso é 


] 
tw 


à, to (o NR e A A E 


H 
ES 
. 


ato 00d 151) 


Desta maneira, O serã selecionado incorretamente como" 
tendo sido a palavra transmitida, quando ê empregado um decodifica- 
dor por distância suave mínima. 

Agora, serã considerado um novo exemplo no qual se assu 
me as palavras código como sendo 0 0 e 1 1, resultando em um côdi- 
go com d=2 e R=0,5. É admitido que a palavra 1 1 é transmitida , 
mas 000111 é recebida na saída do demodulador, ou seja, o 
ruído afetou o primeiro dígito da palavra código. 

As duas saídas de plena confiança são agora , 


0a 1 QT cede trt, de-modor que 


10 


a (0 Odd MI AI) 


aç to 0 ro a ET) 


H 
= 
. 


Portanto, 1 1 serã selecionada corretamente, e um úni- 
co erro foi corrigido, embora o côdigo possua d=2. 

Este tipo de deteção com decisão suave é particularmen- 
te aplicável à decodificação por distância mínima (máximo de veros= 
similhança) para códigos de bloco, e também para códigos convolucio 


nais utilizando o algoritmo de Viterbi ou decodificação sequencial. 


e 


2.2 - ALGUNS PROCEDIMENTOS SUBÓTIMOS 


Nesta seção são descritos algoritmos de decodificação" 
de códigos de bloco que empregam decisão suave e que são subôtimos, 
no sentido que o ganho teórico máximo possível não é atingido. Em ! 
alguns casos isto se deve ao fato do uso de um número finito de re- 
giões de quantização, de maneira a viabilizar a implementação práti 
ca. 

Tais procedimentos, embora subôtimos, representam uma 


melhoria considerável no desempenho com relação a sistemas que uti- 


lizam decisão abrupta. 


2.2.1 - Decisão Suave na Decodificação de Códigos de um Ônico 
Dígito de Paridade 


O procedimento descrito a seguir é devido a Farrell e 
Kalligeros |12|, onde se utilizou como exemplo um código de bloco"! 
(8,7,2) e um detetor com decisão suave quantizado em 8 regiões, cu- 
jos níveis são +0,25V , +0,5V, +0,75V para sinais de +l1V. 
O decodificador deve operar de acordo com a seguinte ' 
sequência: 
1) - avaliar os dígitos de decisão kjk,--«kçk-o, e seus correspon 
dentes "bytes" de confiabilidade. 
2) - Recalcular o dígito de paridade q a partir dos dígitos de 
informação e compará-lo com o valor recebido cj. 
a) Se cy estiver correto, é assumido que não ocorreram erros, 
e a palavra é liberada. 
b) Em caso contrário, inverter o dígito de mais baixa confiabi 
lidade, de modo a corrigir o erro isolado mais provável, e 
liberar então a palavra. 


A curva de desempenho, obtida por simulação em computa 


. 


Ea 


dor, é mostrada a seguir na figura 2.4 abaixo. 

Para relações sinal/ruido maiores que -2,5dB observa-se 
um ganho em relação ao sistema sem codificação, o qual é superior ' 
àquele obtido por um código de Hamming de mesmo comprimento e/ou '! 


eficiência. 


1 - É SNR(db) 
s 
Mo ga! 
DECISÃO 
ABRUPTA 
DECISÃO 
10 SUAVE 


Figura 2.4 - CURVAS DE DESEMPENHO DO ALGORÍTMO DE FARRELL 
E KALLIGEROS. 


. 
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2.2.2 - Decodificação por Distância Suave Mínima 


Decodificação por distância mínima usando decisão suave 


foi empregada em um sistema estudado por Farrell e Munday (1976) ! 


[12 


. É basicamente um sistema de comunicação em HF sem canal de 


retorno. As características do sistema são: 


e) 
2) 


3) 


4) 


3) 
6) 


7) 


baixa taxa de dados, aproximadamente 50 bits/seg. 

Um código produto |36| é empregado, com n = 15x15 = 225 e 
d = 3x3 = 9 codificado e decodificado em cascata. 
Intercalamento de bits para combater a mistura de erros aleatô 
rios e "burst" do canal de HF. 

"Sequence inversion keying" [11] de aproximadamente 100 bits / 
seg, para espalhar o espectro num canal de HF com 3kHz. Para 
tal, uma sequência-m de 15 dígitos foi empregada. 

Modulação ASK, com deteção coerente empregando PLL. 

Correlação dos dígitos da sequência m para determinar se a 
sequência estã invertida ou não (dígito de decisão) eo grau 
de correlação na presença de erros do canal (dígitos de confi- 
ança); o que é equivalente a um esquema de demodulação por de- 
cisão suave com Q=4. 

Decodificação por distância de decisão suave mínima das linhas 
e colunas componentes do côdigo produto, por seleção alternada 
adaptativa das linhas e colunas. 


= 


O código empregado & um código produto, como mostrado ! 


na Ligura 245. 


ne 3 fica 


1 x 11 Dígitos de 


Dígitos de informação |!Ix4| verificação 


de linhas 


pe e, 

— - Verificaçao 
Digitos de verificaçao de sobre 

colunas verificação 


Figura 2.5 - CÓDIGO PRODUTO. 


O desempenho do sistema foi muito bom em testes de rui- 


do. 

Com simulação digital de erros, uma taxa de erros de 
2x10 + nos dígitos da sequência-m foi reduzida para 4x10"* nos 
dígitos de decisão na saída do correlador, e ainda reduzida para 


cerca de 4x1077 depois da decodificação por decisão suave. 
A figura 2.6 mostra as curvas para testes com ruído ! 


branco. 


O ganho do código em ha a aproximadamente f 


12.3dB, enquanto que o ganho teórico máximo & cerca de 14.7dB. A di 
ferença não é grande tendo em vista que esta degradação possibilita 


a realização prática. 


2.2.3 - O Algoritmo de Harrison 


Pode ser mostrado |40 | que o processo de decisão abrup- 
ta num canal gaussiano resulta numa degradação no desempenho de 
aproximadamente 2dB comparado com o receptor Ótimo. Segue-se que |, 
para este canal, o máximo de melhoria a ser esperada por introduzir 


decisão suave no sistema, é equivalente a esta degradação. Harrison 
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|17| propôs um algoritmo de decodificação empregando decisão suave" 
que permite um ganho de aproximadamente 1dB, o que é um bom resulta 
do tendo em vista a relativa simplicidade de implementação do mesmo. 
Neste esquema, a palavra recebida consiste da informação básica da 
decisão abrupta com "marcas" de confiabilidade (indicando quais as 
prováveis fontes de erro). Esta informação de confiabilidade dispo- 
nível ao decodificador indica simplesmente quais os bits na palavra 
código que são os mais prováveis de terem sido afetados por erros. 

Utilizou-se um código de Hamming (7,4,3) como exemplo ! 
para aplicação desta técnica. Desde que a complexidade do decodifi- 
cador cresce rapidamente com o aumento do número de níveis de quan- 


= 


tização, um esquema com 4 níveis ê descrito, 


A 
A/2 
H 
DO Cnnrapeçdia 6:46 foi DUM gire UR palc (6) (079 H > Alta confia- 
L bilidade 
q L > Baixa confia 
bilidade 
-=J 
-A/2 


Figura 2.7 - REGIÕES DE CONFIABILIDADE. 


A posição dos limiares relativos a confiabilidade é mui 
to importante nas considerações de informação de confiabilidade. Se 
a tensão J é muito pequena, poucos bits recebidos serão de baixa! 
confiabilidade, transportando pouca informação ao decodificador. Em 
contrapartida, se JU é uma tensão elevada, praticamente todos os 
bits recebidos serão de igual (baixa) confiabilidade, o que é equi- 
valente a decisão abrupta. 

Existe uma posição Ótima, para cada código corretor, a 


qual maximiza a informação de confiabilidade entregue ao decodifica 


A 


ya 


dor. As curvas de desempenho obtidas para 3 limiares de decisão di- 

ferentes são apresentadas na figura 2.8, onde ê observada a exis - 

tência deste valor ótimo para J entre 40% e 50% da amplitude dos 
pulsos binários transmitidos. 

O procedimento proposto consiste basicamente dos seguin 
tes passos: 

1) - Se a palavra recebida possui 0,1 ou mais que 2 dígitos com bai 
xa confiabilidade, o decodificador atua desprezando a informa- 
ção probabilística, ie, empregando decisão abrupta. 

2) - Se a palavra recebida tem 2 dígitos com baixa confiabilidade , 
então a sindrome s é calculada: 

a) Se s = O, então a palavra e suposta correta e ê liberada. 
b) Se s indica um erro em uma posição de baixa confiabilida- 
de, então a mesma é corrigida, e a palavra liberada. 
c) Se s indica um erro em uma posição de alta confiabilidade, 
então dois dígitos de baixa confiabilidade são invertidos e 
a sindrome é recalculada. Se o novo valor é s = 0, é consi 
derado que dois erros foram corrigidos, e em caso contrário, 
ocorreu um erro numa posição de alta confiabilidade e o der 
codificador deve operar normalmente. 
Desta maneira o desempenho do decodificador com decisão 
suave deve ser pelo menos tão bom quanto o receptor equivalente em- 
pregando decisão abrupta, para todos os valores de relação sinal / 


ruído. 


2.3 - DECODIFICAÇÃO POR DECISÃO SUAVE PARA CÓDIGOS DE BLOCO 


USANDO UMA TRELIÇA 


Serã descrito um algoritmo introduzido em 1978 por '! 


Jack Wolf [38], o qual faz uso de decisão suave na decodificação por 
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Figura 2.8 -— ALGORÍTMO DE HARRISON - CURVAS DE 
DESEMPENHO. 


Ai 


máxima verossimilhança de qualquer côdigo de bloco linear (n,k,d) so 
bre GF(q), e que pode ser implementado usando uma treliça tendo não 


mais que Pi estados. 


2.3.1 - Considerações Gerais 


O algoritmo apresentado faz uso da técnica de decisão ! 
suave, O que estã associado ao fato de que emprega números reais (eg, 
a saída analógica de filtros casados para os sinais) relacionados ! 
com cada símbolo correspondente da palavra código. 

Serã mostrado que a decodificação de côdigos de bloco ! 
usando uma treliça é de grande interesse na decodificação de códi - 
gos com alta eficiência, visto que a complexidade da treliça é limi 
tada superiormente por uma função do número de símbolos de paridade. 
Também, algumas simplificações resultam quando são considerados cô- 


digos lineares particulares, tais como códigos cíclicos ou códigos" 


produto |38|]. 


2.3.2 - Treliça Associada a um Código de Bloco Linear 


Seja a matriz de verificação de paridade [H] de um cô- 


digo de bloco linear (n,k,d) em GF(q) considerada sob a forma 


[H] = E h,.--h,|, onde cada vetor coluna h, possui (n-k) 'co- 
a = E 

ordenadas em GF(qg), ou seja, [8 |- [as hecho 4) d=l,2,...,n. 

Considere um diagrama de ârvore como uma coleção de 


nôs (estados) interligados entre si por ramos unidirecionais. Cada 
estado a uma dada "profundidade" k & indicado por uma (n-k) -upla 
S; (k) com elementos em GF(g). As linhas são traçadas da profundida 
de k para a profundidade k+1, considerando os nôs (ts, (k6)) a uma 


dada profundidade gerados como indicado no diagrama da ELlgura 2.9, 


RÁ fm 


(Profundidade) K=0 K=1 =2 =3 


So+ã h1l 


. 
. 
. 
. 

o 


Ss 


+qh1+qh2 


Figura 2.9 - DIAGRAMA DE ÁRVORE PARA UM CÓDIGO 
DE BLOCO. 


Deste modo é possível verificar que a uma dada profundi 
dade k existem os seguintes nôs: 
Ea 
s RR) 01 
Sq Das, = & <k 
SO as dh, Das, + 


Ix 
Ês 
A 


IA 
+ 
A 
[25] 
A 
x 


etc 
De uma maneira resumida, é possível escrever que a uma 


dada profundidade k existem os nós: 


k 
Sp (k) = So (0) O; h;. com a; e Gr(g) (2-1) 


Cada estado ê identificado por uma (n-k)-upla em GF(g), 


de modo que o número de estados em qualquer profundidade não pode ! 


exceder q , O número de  (n-k)-uplas distintas com elementos em 
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GF (q). Deste fato segue-se que & possível fazer uma modificação na 
figura 2.9, que consiste em juntar os nôs indicados pela mesma g 


(n-k) -upla, resultando em um diagrama de treliça. 


Definição 2.3.2 - Uma treliça associada a um código de bloco em 
GF(q) & uma coleção particular de nôs (ou estados) interligados '! 
por ramos unidirecionais, e agrupados em conjuntos indexados por 
um parâmetros k, k=0,1,2,...4,n. 

Um nô indexado por um dado valor de k & dito estar 
a uma profundidade k. As linhas serão traçadas entre certos pares" 
de nôs a uma profundidade k e a uma profundidade k+l, para ' 
k = 0,1,2,... n=-l; com a direção do ramo da profundidade k para a 
profundidade Kk+l. Para qualquer profundidade Xk existirão no mãxi 


je 


mo q” nós, os quais serão identificados por (n-k)-uplas, S; (X) 


com elementos em GF(q) para certos valores de Ji. Todas as Pica ! 
(n-k) -uplas são ordenadas de 0 a E als TR com O referindo-se a 
(n-k) -upla toda zero. S, (k) é interpretada como a i-êsima 4 
(n-k) -upla desta lista ordenada. Desde que nem todas as (n-k) -uplas 
podem corresponder aos nôs a uma profundidade k, deve ser considera 


do um conjunto Pp, subconjunto dos inteiros E e AR ê 
correspondendo a essas (n-k)-uplas que estão relacionadas a nós a 


uma profundidade XX. 


2.3.3 - Construção da Treliça 


A treliça associada a um código possui n+l profun- 
didades, desde k=0 a k=n, e possui a cada profundidade E esta 
dos possíveis. Cada nó ê caracterizado pela profundidade k e pelo 
estado A ER RD ia de modo que cada nô grigina q ou- 
tros. Cada nô em qualquer profundidade k & obtido como combinação 
linear dos vetores h, correspondentes a profundidades menores que 


L 


A 
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k, de acordo com a equação (2-1). É importante notar que S, (x) não 
& univocamente definido, pois diferentes caminhos (no máximo q) po- 
dem levar a um mesmo estado. 

O número de nôs que são definidos a uma profundidade Xk 
& dado por a, onde LI é o número de vetores coluna (né li - 
nearmente independentes. 

A construção da treliça associada a um código de bloco! 
& feita de acordo com o procedimento descrito a seguir: 
1) - Na profundidade k=0, a treliça contêm somente um nô, chamado! 

So (0), a (n-k)-upla toda nula. 

2) - Para cada k = 0,1,...,n-1l, a coleção de nós na profundidade ' 
k+1 & obtida a partir da coleção de nôs na profundidade k pe 
la fôrmula: 


Sp(k+1) = S,(k) Da, hi, iePyr apeGr(g) (2-2) 


Para cada icePpo, as linhas de conexão são traçadas entre o nô 
8; (k) e os q nôs na profundidade k+l, gerados empregando a 
equação (2-2) acima. Cada linha é determinada pelo valor parti 
cular de a, que gerou Sp(k+l) a partir de S; Kk). 

3) - Expurgar quaisquer nôs que não tenham um caminho para o estado 
todo zero na profundidade k=n, e retirar todas as linhas tra- 
çadas para estes nôs expurgados. 

Hã uma correspondência biunívoca entre cada palavra cô- 
digo e a sequência de a 's, sobre qualquer caminho, do nô todo ze 
ro na profundidade q para o todo zero na profundidade mn. Existem ' 
qÉ caminhos distintos atravês da treliça, e cada um deles corres - 
ponde a uma palavra código. Assim, a treliça fornece um mêtodo com- 
pacto de catalogar todas as palavras código; cada palavra código . ! 


distinta é representada por um caminho distinto na treliça. 


- 


Exemplo 2.1 - O procedimento descrito & melhor ilustrado neste exem 


plo, onde é considerado o côdigo binário (5,3,2) introduzido no o 


=Ãio 


exemplo 1.1, cuja matriz de verificação de paridade foi determinada 


no exemplo 1.2 como sendo 


[8] To JQREL O 
H] = a) = |h; Bb; ho b, À 
ca Aqua pi Oi 


Neste caso, n=5, k=3 e q=2, de modo que q Es e 
Pp ElO0;1,263%. 
Os nôs são associados aos números inteiros na forma in- 
dicada: 
0 0 1 1 
> 0 ag d > 2 e + 3 
(o) PR PM e 1 


A treliça para este código ê apresentada na figura 2.11. 


Prof. 
Ratados Ke 


0 K=1 K=2 K=3 K=4 K=5 
(8) 
(MD . 
(D + 
O 


Figura 2.10 - TRELIÇA NÃO EXPURGADA PARA O 
CÓDIGO (5,3,2). 


NR cin K=0 100 Ke2 K=3  K4 Kes 
(9) 
(De 
(D+ 
(De 


Figura 2.11 - TRELIÇA PARA O CÓDIGO DE BLOCO (5,3,2) 
(EXPURGADA) . 


sab 

A referência |25| ES pe algoritmos, "oscquais 
permitem uma redução na complexidade da busca em uma treliça asso — 
ciada a um código de bloco. Entretanto, ê evidente que deve existir" 
um compromisso entre a complexidade e eficiência dos sistemas que 
empregam tais algoritmos, o qual ê analisado atravês de simulação , 
na referência citada. Alguns destes procedimentos serão apresenta - 


dos na seção 2.6. 


2.3.4 - Decodificador de Viterbi usando a Treliça 


A decodificação & executada baseada na ênupla recebida! 
£ + com componentes reais. É assumido que não hã interferência inter 
simbólica, de modo que a j-ésima componente de r depende somente 
da j-êsima componente da palavra transmitida c. 

Também é assumido que a contribuição do ruído em cada 
uma destas componentes ê descrita por variáveis aleatórias estatis- 
ticamente independentes N;r com fdp Ea ( jd = resido” Então 
o logarítmo da função de verossimilhança, dado a palavra código ' 


transmitida, &ê da forma 


n 
log £ lo (rle) = os log (eu, “estoy | = 
n 
- E, Bs(reep) = Ste) (2-3) 


Para uma dada sequência recebida r, & óbvio que uma 
decodificação com máxima probabilidade a posteriori é aquela que 
encontra a palavra código c que maximiza Z(c). Seria emprego da 

é ' k 
força bruta verificar por exaustao, ie, tentando todas as q pos 
hero Pe 
siveis palavras côdigo. 
A seguir é descrito um algoritmo recursivo pelo qual ! 


muitas palavras côdigo podem ser descartadas no processo de decodi 
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ficação quando na determinação da palavra código que maximiza Z(c). 


A cada nô S; (K) na profundidade k, deve ser associado 


um número real VíS; ()) de acordo com as seguintes regras: 
1) - Para a profundidade k=0, faça V(Sg (0)) = 0. 
2) - vL € Pr” forme V(S, (k+1)) a partir de VíS; ()) da seguinte 
maneira: 
V(Sp (k+1)) = Max (US; +24 Msg) (2-4) 
a ,EGF (q) J 


com iePp! onde E ao P & um subconjunto de Índices tal 


que, para algum “4 de GF(q) tem-se 


Spk+1) = s (6) O “ ba (2-5) 

3) - Retenha somente aquele caminho para o quala fórmula (2-4) forne 
ça o máximo. 

4) - Em k=n, a sequência de a;'s obtida simplesmente ligando o ca- 


minho do estado todo zero na profundidade k=0 até o estado to 
do zero na profundidade k=n, corresponde à palavra código Te) 


que maximiza Z(c). 


2.4 - DECISÃO SUAVE NA DECODIFICAÇÃO QUANDO A FONTE TEM DISTRIBUIÇÃO 


DE PROBABILIDADE DESCONHECIDA 


Serã apresentado um procedimento que faz uso de coefici- 
entes probabilísticos condicionais de confiança como medida de con - 
fiabilidade, resultando em uma regra de decodificação de máxima ve - 


rossimilhança condicional. 
2.4.1 - Noções Preliminares 


Considere um sistema de comunicação no qual os símbolos" 
da fonte binária tem uma distribuição de probabilidade desconhecida! 


e são transmitidos em blocos de n dígitos atravês de um canal sem 


. 
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memória, caracterizado por uma função de verossimilhança f(y, |j),'! 
onde E LOG; y; € RE 

Uma vez a palavra binária transmitida e uma sequência '! 
x = (yy'-::+y,) recebida, o uso do método de Neymann-Pearson clãs- 
sico |13] para cada y, com um nível de confiabilidade prescrito , 
resulta em um dígito binário x,- À decisão x; não considera qual 
quer medida de conclusividade relativa a uma observação particular"! 
, 


y.: pois diferentes Y; 


á s podem ser decodificados em um mesmo x,, 


a é 
com as mesmas probabilidades de erro tipo I e tipo II, podendo cada 
uma das decisões possuir diferentes graus de conclusividade. 

A maioria dos procedimentos de decodificação que empre- 
gam decisão suave, tratam de casos quando a distribuição da fonte é 
conhecida, de modo que & possível ser usada a probabilidade a poste 
riori como medida de confiabilidade. Contudo, & conhecido que em 
vários problemas práticos esse não é o caso. Wolfenson e Rocha |39| 
propuseram fazer uso de coeficientes de confiança condicional de 


Kiefer |21|, os quais têm interpretação frequentista, como uma medi 


da de confiabilidade. 


2.4.2 = Decodificação usando Confiança Condicional 


O uso do lema fundamental de Neymann-Pearson permite de 
codificar com probabilidades de erro prescritas, cada componente da 
palavra recebida y, resultando num vetor x. Como discutido, este 
procedimento não expressa nenhuma conclusividade levando em conta ' 
uma observação Yi particular. 

Denotando por ae 8 as probabilidades de erro tipo 
e II, a conclusão provida pelo método é simplesmente que cada Yi 
deve ser decodificado em x, com probabilidades de erro a e 8 . 


Kiefer |21| propôs um esquema de atribuir uma confiabilidade a cada 


decisão x; em termos de um coeficiente de confiança condicional. 


. 


ao ia 


Considere uma partição Il do espaço de observações JR, 
onde I = tus) para je (0,1), se S (S um conjunto de Índices). 


É também definido IÉ através da relação nê=n5 Ur . 
â , 


Observando um dado Y; € R, é procurado o elemento 


1 e II no qual Yi; caiu, e assim Y; é decodificado em 5 com ! 


* 
confiança condicional dada por Rj, = P (dec |1ºº,5). 


É assumido ainda que W & escolhida de tal maneira que 
0 0 
pídec O |11º',0) = Pídec 1/19 ,1) (2-6) 


Rio sendo uma probabilidade, pode ser interpretada em termos de 


frequência relativa. Repetindo o experimento independentemente um 


número muito grande de vezes de modo que se tenha um grande número! 


: so ; e 
de ocorrências de IT |, Roo fornece aproximadamente a fração das 


so E ' RR 
vezes que IT ocorreu e nas quais foi tomada uma decisao correta. 
Isto fornece, portanto, uma maneira objetiva de atribuir uma medida 


de confiabilidade a cada decisão NP. Assumindo que Y= (yyrecceYh) 


foi recebido e que seja verificado y; E qe para cada 5 


i =1,2,...,n, Segue-se que cada pai pode ser decodificado pelo ' 


procedimento convencional NP com ae 8, e a medida de confiabilida 


de da decisão é dada por Rs;º À questão agora ê como combinar as 


confiabilidades Roi! i=1,2,...,n de modo a atribuir uma medida" 


de confiabilidade à sequência decodificada x = (xy rece rXdo Se x 


é a sequência decodificada, levando em consideração que o canal ê 


sem memória, ê possível escrever 


n 


si 
= sa 


À n 
pídec x| [11º*] c x)=T Pídecx, |X', 


n 
=T Ro = R(g) (2-7) 


* Veja lista de símbolos e abreviaturas. 
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Desta maneira & associada uma confiabilidade dada por 


n É 
R(x) = T Ras à sequência inteira decodificada. 


2.4.3 - Exemplo: Aplicação a um Canal com Ruído Branco 


Gaussiano 


O método de decodificação apresentado é uma solução de 
compromisso entre minimizar a probabilidade de erro por dígito e es 
colher a palavra código mais provável. É possível visualizá-lo como 
um procedimento em dois passos, onde primeiramente o critério NP é 
usado, com probabilidades de erro prescritas, e então a estrutura ' 
do código & empregada juntamente com a medida de confiabilidade, de 
modo a encontrar a palavra código mais provável, condicionada E] 
observação. 

Nesta aplicação será considerado um canal perturbado 
por ruído branco gaussiano com média | nula e variância unitária. A 
fonte bipolar transmite +1V ou -1v correspondendo aos símbolos 1 
e 0, respectivamente. 

Também & assumido iguais probabilidades de erros oq =8 


a decisão NP é y;>0 com e = BI= Ot) = 0,16. 
0 


A partição Il é escolhida de modo a permitir a máxima 
variabilidade nos coeficientes de confiança. No exemplo, S & esco- 
; , : : s ? 
lhido como o conjunto dos possíveis valores de Rçr ie, I será in 


dexado de modo que: R$ = s.- Portanto, 


Pídec 3ly,3) = Pídec j,y/j) (2-8) 


Ply|3) 


bs) 
H 
u 
l 


ou seja, 


R=se= (2-9) 


s 2 A 
é Pídec 1-j,y/5) 


Pídec 5,y]5) 


Assim, s = : e 
T ue NCDA 
; 1 + er y 
nº = (£ 1/2 


Suponha, por exemplo, que 


de forma que o procedimento 
com coeficientes condicionais 


lidade R(x) = 0,207. 


Observe que devido a escolha de 


codificado em 1 ou O 


E es 


1º |, & escolhida como 


tn (s 1.3) para 80 1/2% (2-10) 


y = (-0.2,-0077-0:7,0,29121) 
NP bit a bit resulta em x=(0,0,0,1,1), 


(0.6,0.8,0.8,0.6,0.9) e com confiabi 


Il, +0.2V pode ser de- 


mente. Supondo que 


no contexto acima, 


o 


ê 


com medida de confiança 0.6 ou 0.4, respectiva 
código (5,3,2) descrito no exemplo 1.1 é usado 


possivel fazer uso do diagrama de treliça asso 


ciado (veja exemplo 2.1, figura 2.11) para encontrar a palavra códi 


go de maior confiabilidade. 


Os passos da decodificação estão esquematizados abaixo. 


00 0.6 0.48 0.384 
01 0.08 0.096 
10 0.12 0.096 
11 -4 0.32 0.256 
(a) Cb) (c) 
0.057 0.138 
0.153 
(d) (e) 
Figura 2.12 - PASSOS NA DECODIFICAÇÃO DO CÓDIGO (5,3,2) 
USANDO MÉTODO DE WOLFENSON-ROCHA. 
Observe que um erro foi corrigido no primeiro dígito 
(R = 0.6), pois, uma decisão abrupta estima a palavra transmitida 
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como (0,0,0,1,1) com R(x) = 0.207, citada que a saida do decodi- 
ficador estima (1,0,0,1,1) com R(x) = 0.138. Deve também ser obser 
vado que a posição alterada foi a de menor confiabilidade, e que 
uma mudança na quarta posição (R$4 = 0.6) não é permitida pela es 


trutura do código. 
2.5 - DISTÂNCIA GENERALIZADA 


Serã mostrado que distância generalizada constitui uma 
extensão da idéia de apagamento discutida anteriormente, a qual per 
mite uma maior flexibilidade no uso da informação de confiabilidade. 
O emprego de distância generalizada permite que algorítmos algêébri- 
cos de decodificação possam ser prontamente adaptados para fazerem" 
uso da decisão suave. Desta forma, algumas características de deco- 
dificação probabilística são introduzidas sem sacrificar a atrativi 
dade dos procedimentos algébricos. 

Serão analisados dois algoritmos para decodificar por 
distância mínima generalizada, sendo o segundo procedimento descri- 
to, proposto pelo autor. Para o primeiro procedimento, introduzido! 
por Forney |14|] em 1966, foram tambêm desenvolvidas -cotas exponenci 


ais para a probabilidade de não decodificar corretamente. 
2.5.1 - Medidas de Distância 


Os algorítmos de decodificação são formulados, geralmen 
te, em termos da escolha da palavra-côdigo mais próxima da palavra! 
recebida, onde a "proximidade" & definida por alguma medida de dis- 
tância. Duas importantes medidas de distância são apresentadas nes- 
' 


ta seção. A distância de Hamming, já introduzida na definição 


1.2.2bp serã' reapresentada deruma maneira mais formal. 
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Definição 2.5.1a: - A distância de Hamming Dy(£,g) entre duas pala 
vras fe g é simplesmente o número de posições em que elas diferem; 


n 
ca Ag (f; 194) 


ley D. (Ed) 
nes i=1 


Onde 


A (Es sg) = 
Egas se E, Í g 


Suponha que o receptor decide, para cada amostra recebida, 
qual das q letras foi a mais provavelmente enviada, uma escolha ! 
que corresponde a um elemento do GF (q) rotulado por re Esta esco - 
lha é então associada a uma classe de confiabilidade e; de acordo! 
com a certeza que se tem da mesma ter sido correta. A cada classe es 
são associados-dois parâmetros Bo3 e Bai tais que 0585958095) a 
Contudo,serã mostrado que somente a quantidade as Sab qm Bs Jena 
mada de peso da classe e; ê significante. Vê-se claramente que 


(vj) & verificada a relação 0 <a, 


hos 
jc 


Definição 2.5.1b - Distância generalizada Dç tz. £) entre a palavra! 


recebida r e a palavra código f & definida como sendo 


A 
Da(x,£) o E dg (r E Epa 
i=l 
Onde 
Go SS E, SH, , Ty ER: 
e E de 5 
Aalrjrf,) = 

Ba) se r, Ed É, XL [5 e 


Deve ser observado que esta distância não constitui uma 
métrica |2 |, enquanto que Dy representa uma métrica. 

Em geral, as classes para as quais a=l serão referidas 
como classes de plena confiança, e aquelas para as quais a=0 como! 
classes não confiáveis. Valores intermediários de A correspondem! 


a níveis intermediários de confiabilidade. Também serão ordenados os 


números das classes em ordem decrescente com relação ao seu peso, de 
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maneira que se 5 <k, então as 2 


2.5.2 - Decodificação por Distância Mínima Generalizada I 


Algumas das propriedades da distância generalizada se - 
rao discutidas a seguir, as quais permitem o estabelecimento de um 
decodificador usando distância mínima generalizada. Cotas exponen - 
ciais de erros para o algoritmo estabelecido serão também desenvol- 


vidas. 


2.5.2.1 - Propriedades da Distância Generalizada 


Inicialmente suponha que durante a transmissão de uma 


palavra f de um código de bloco (n,k,d), o número de letras rece- 


- 


bidas corretamente (r; = Ê£,) e postas na classe “3 e o 


Dos 
número das recebidas incorretamente (r; al Ê,) e postas nesta classe 
5 ej 
Teorema 2.1 - Se £ é enviada e Ao; * Das são tais que se verifi 
= + (Il+a. bo rd a RE 

ca : (a dba ( 043) , então Do (x, E) Lg) para 
todas palavras código g f £. 

Prova - Pela definição de distância generalizada segue-se que 


Do tz,£) (2-11) 


H 
mM 
(= 
To 
Q 
“u 
Q 
q. 
To 
m 
qu 
D 
u 
ips 


Para qualquer palavra código gr E serão considerados 
três conjuntos de Índices, de acordo com 
E) se Tt Ij 
Sie : 
cj £ 


se £f. *Z g 


ENT pa 


O número de letras |s | e | | dos conjuntos s. e 
cj eJ cj 


s . Satisfazem as relações |s 


z e |s | nm .« Então: 
| ejl = “ej 


| 4 Reg 
c7| = “e3 


Açlr;19;) > 0 = dy(g;,f;) ies, (2-12) 
As(r;19;) = Bos = d(g;rf,)-1 + Bos LE 8ç3 (2-13) 
dA(r;19;) > Bo3 = Aglg, 1 f,)-1 + Bc3 ie So3 (2-14) 


Operando com as equações (2-12) a (2-14) acima, obtêm - 


se: 


Dç(z.9) > Dy(£rg) — (0-0) [Ses | + (18,5) [Se3 | 
fabio RA 2 (ctg nes + (-859)ne5) (2-15) 


e oO o “)n 


jli ipõ > -B..+ : 
Utilizando a hipótese d o eg Ea ei! Pes 


para enfraquecer a desigualdade (2-15), tem-se 


Dçtz.9) > : (Pos Bia Aga 8.3) = Dç(z:£) 
C.Q.D 
É óbvio que se o critêrio de distância mínima generali- 
zada for utilizado como critêrio de decodificação, não serão cometi 
dos erros nos casos onde Nç3 e ho3 são tais que a desigualdade do 
teorema 2.1 & satisfeita. 
No caso especial em que hã somente uma classe qq7+=L, 
na é o número total t de erros em r, de maneira que o teorema" 
torna-se Dç (xr £) < DçtL rg) vg * £ se 2t <d. Se for permitida a 


inclusão de uma classe de apagamentos «à, = 0, cujo número total é 


s = Nç2 + Do! o teorema pode ser reescrito como 


Dç(r.£) < Dglxeg) Vg É £ se 2t+s<d. 


io 


Estas são condições familiares que asseguram a não ocor 
rência de erros quando for utilizada a decodificação por distância! 
mínima. 

Definição 2.5.2a - É dito que a palavra código estã "dentro da dis- 
tância mínima da palavra recebida" quando a desigualdade do teorema 
2.1 é aparentemente satisfeita. 

O teorema mostra que hã no máximo uma palavra côdigo '! 
dentro da distância mínima de qualquer palavra recebida. Assim, se 
o decodificador encontra alguma palavra código dentro da distância 


mínima da palavra recebida, ele pode imediatamente anunciar que es- 


ta é a palavra decodificada. 


2.5.2.2 - O Algoritmo de Forney 


Forney Jr.|14] propôs um algorítmo prático que usa dis- 
tância mínima generalizada como critério de decodificação, e que 
funciona se hã uma palavra dentro da distância mínima da palavra re 
cebida. Aqui é assumida a existência de algoritmos de decodificação 
que são capazes de manipular com apagamentos e erros, e que operam! 
se o número s de apagamentos e o número t de erros são tais que 
A co O = Radar o UA 

Suponha, para este propósito, que as letras recebidas ' 
em algumas classes são consideradas como apagamentos, e as restan - 
tes como completamente confiáveis. Denotando por 


- 


E = Ule, é não confiável), e R = ES =(tile, é de plena 
confiança) 
é tentado decodificar a palavra recebida pelo algoritmo de erros-e- 


apagamentos, o qual estarã apto para decodificar se 
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+ - 
2 Mia ns e Nos) Sds (2-16) 


Es fics 
jer jeE 

Se isto acontece, verifica-se, então, se a palavra coódi 
go obtida estã dentro da distância mínima da palavra recebida. Caso 
ela esteja, então foi determinada a única palavra dentro da distãân- 
cia mínima. 

Não hã uma única atribuição provisória de apagamentos ! 
para a qual este método tem sucesso. Contudo, o teorema a seguir e 
seu corolário mostram que um pequeno número de tais tentativas apre 
sentam êxito em encontrar a única palavra código dentro da distân - 
cia mínima, se existe uma. Como anteriormente, as classes são consi 
deradas ordenadas de acordo com ordem decrescente de confiabilidade. 

Seja J o número total de classes, e um vetor de JU dá 
mensões a definido como a= (Ay rAgre cond). São também definidos 
os conjuntos R, = (5lj <b) e E, = (9]9 > b+r1Y com Q< bx dg. 

É considerado que Sr um vetor J-dimensional com 1's 
nas primeiras b posições e zeros nas demais, representa uma atri- 
buição provisória correspondente a R = R, eE=E- A idéia do teo 
rema 2.2 abaixo é que a estã dentro de uma envoltôria convexa, cujos 


pontos extremos são a, 's, enquanto a expressão f(a) a ser defini 


b 
da é uma função linear de a, que deve assumir valores mínimos em 


algum ponto extremo, ie, para uma das atribuições provisórias 4 


J 

Teorema 2.2 - sê Fios gos q + (ta j)no;h a “4 2 4 para 
j < k, existe algum inteiro b tal que a desigualdade : 

b J 
2 £- nn.c+rk (n .+n..)]) < d. é satisfeita. 

3=21 2 s=be 01) E 
Prova - Inicialmente é definida uma função do vetor a, 

[oo ul ” 
f(a) 2 E j(1-0)n.. + (lta)n ) < (2-=17) 
(a) as: do E) a] 
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a qual satisfaz as seguintes propriedades: 


i - £ & uma função linear do vetor J-dimensional a 
b J 
ii= note que f(a,) =22Z n.+2Z Gra tm) 
> j=1  s5=p1. SJ e) 
O teorema é provado supondo que (vb) E(g) >d, 
0 <b<J, e exibindo uma contradição. 


Para isto,sejam An = Ll=-a 


0 1 
Ab “0 — 4 | ese ts E 
As “As 
É facilmente verificado que 0 S A < 1 para q 
J 
& Dilma pd “6 que z A, =1. 
b=0 
J 
Ora, a=kL A a, demodogqueéê possívelescre- 
ppp ESB 
ver 
J J J 
Eta) = ECZoA ls CÊ AM EM) SAL A Sd 
io bao POD qeg bh b=o P 
CONTRADIÇÃO ! 
Portanto, conclui-se que E (a) < d para pelo me - 
nos um valor de b, 0 <b<uJ. €.0.D. 


Este teorema prova que se hã alguma palavra código! 
dentro da distância mínima da palavra recebida, então deve haver al 
guma atribuição provisória, a qual permite que o decodificador de 
erros-e-apagamentos tenha sucesso em decodificar de acordo com o 


critêrio de distância mínima generalizada. Mas um decodificador de 


erros-e-apagamentos tem êxito se e somente se existem aparentemente: 


Nenhum erro e d-l apagamentos 
Um erro e dos apagamentos 
=| d=1 
to erros e do = 2tg-1 apagamentos, onde to me 3 . 
Corolário - Bastam totIls< E tentativas de possíveis atribui 


ções provisórias para que se tenha sucesso em decodificar qualquer" 


palavra recebida que estã dentro da distância mínima pelo critério" 


Eras (E 


de distância minima generalizada, não importando quantas classes de 
confiabilidade existam. 
Segue-se, portanto, que o número máximo de tais proces- 


sos & somente proporcional a d. Ademais, muitos dos processos (tal 


vez todos) podem ter êxito, de modo que o número médio de proces - 


sos pode ser apreciavelmentemenor que o mãximo. 


2.5.2.3 - Cotas (exponenciais) de Erros para o Algoritmo de 
Forney 


Para o algorítmo proposto por Forney |14]| existem desen 
volvidas cotas "apertadas" para a probabilidade de não decodificar" 
corretamente, onde "não decodificar corretamente" significa decodi- 
ficação incorreta ou incapacidade de decodificar. Este último even- 
to ocorre quando não hã palavra código dentro da distância mínima ! 
da palavra recebida. 


Considere uma variável aleatória yjra qual para cada 


letra transmitida assume valores: 


1-6. se a letra é recebida corretamente e posta em es 
à pe 
be se a letra ê recebida incorretamente e posta em cs 


Assumindo um canal sem memória, estas variáveis aleatô- 
n 
rias ty;h resultam iid. 


Sejam PojPrty;=(10)) Je Poj Prtyj=(l+a,)) (2-18) 


A função geradora de momentos |8g | g(s) para a variã- 
vel aleatória y, & dada por 
1-a.) s(l+a.) 
RD j 


Des 
g(s) = Efe ?*) = E Poj € +p.e (2-19) 
5 
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e a função geradora semi-invariante |8 | correspondente, expressa ' 


por u(s) = êng(s). 


Jã é sabido que não haverã erros na decodificação se 
Nç5 e No5 são tais que a desigualdade do teorema 2.1 é satisfeita, 
n 
o que ê o mesmo que requerer que bx 6) < d. Deste fato segue-se ' 
i=1 
que 
n 
pr(nde) = prí E y; 2 dj (2-20). 
i=1 


A cota de Chernoff | 8 | é uma cota exponencialmente , 
apertada para a probabilidade de que a soma de variáveis aleatórias 


iid exceda um certo número; no caso, a cota resulta em 
pr(ndc) < exp -n |s6 - u(s)| (2-21) 


válida vs > 0, onde S =d/n. 


É possível tornar a cota mais apertada maximizando o ex 


poente pela escolha de s e dos a.'s: 


a] 


E(6) = max |s6 - u(s) | (2-22) 


S/a, 


Primeiro, a maximização serã feita sobre Aê e desde 


que u(s) = tng(s), isto & realizado minimizando g(s). 
3 E (Lg) s(ltas) 
a (a) 8 poe T+ Pose =0 (2-23) 
J 
os a's ótimos são então determinados a partir de 
(2-23) acima, lembrando da restrição O < as < 1, o que resulta em 
0 se E 
j é 
a L;/25 se Q < b; <2s (2-24) 
1 se La a As 
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CR prír, Correta] r. e c.) 
onde Li; = [te RS a SR 
) Pe3 prír, incorreta|r, E es) 


Assim, a atribuição Ótima dos pesos envolve apagamento" 
de qualquer recepção para a qual a probabilidade da escolha correta 
p seja menor que 1/2 ; considera plenamente confiável qualquer re- 
cepção para a qual p seja maior que um limiar (que depende de SS), 
e para valores intermediários assume um valor proporcional ao Log 
da razão de verossimilhança tn p/(1-p). Portanto, exceto pelas li- 
mitações nos extremos, a decodificação por distância mínima genera- 
lizada é um mêtodo que utiliza o Log da razão de verossimilhança em 
esquemas de decodificação algébrica. 

Suponha definidas as classes “JeR se L. > 2s 

jeE se L. <o0O 
jeG noutros casos 


Segue-se, então, que 


ÉS) TE 2s s s = 
Iopt e Po (s) + e Py(s) + e pç (8) + Po(s) (2-25) 
onde, 

Pos) = E po; + Pçls) = E Pej (2=26) 
jeR E jeR 

(8) = 2 ip. + o) 2-27 
Px dt Pas * Pá ( ) 

= E 2Yp-a R 2-2 
Pgts) eU Po3 Po5 ( 8) 
Finalmente, assumindo Hopt (8) = tn Iopt (8) e tendo em 

vista a equação (2-22), tem-se que 

E(6) = max |sô - VUopt 8) | (2-29) 


s 


Agora otimizando pela escolha de s, observa-se que 


aGã= 


para o s ótimo verifica-se a relação 


O mit (2-30) 


= 
Como Vopt 8) é monotonicamente crescente com s, pois, 
u(s) & convexa | 8|, a equação (2-30) acima admite solução com 


s > 0 se e somente se 6 > (0). É concluido, portanto, que a 


,' 
Hopt 
menor distância mínima (de Hamming) para a qual a decodificação por 
distância mínima generalizada pode operar ê dada por niop+ (OU - 

2.5.3 - Decodificação por Distância Mínima Generalizada II 


Neste seção são feitas algumas considerações adicionais 
sobre distância generalizada, apresentando um algoritmo de decodifi 
cação por distância mínima para códigos de bloco, o qual faz uso da 
treliça associada ao código. 

O algoritmo aqui apresentado estabelece um procedimento 
equivalente ao proposto na seção anterior, no sentido de que também! 
escolhe a palavra código mais initaião da palavra recebida, quando ê& 
utilizada como medida de proximidade a distância generalizada. Deve 
ser mencionado, entretanto, que existem algumas dificuldades práti- 
cas na implementação do decodificador descrito pelo teorema 2.2 da 
seção anterior. As tot 1 tentativas de decodificação do corolário! 
deste teorema podem resultar em até tot 1 palavras código diferentes, 
A classe de algoritmos apresentados por Forney sugere o uso de téc- 
nicas algébricas para gerar um certo número de palavras código prô- 
ximas em algum sentido da palavra recebida. O problema crucial fé 
encontrar uma técnica eficiente de gerar o Conjunto de palavras có- 
digo, o qual tenha alta probabilidade de conter a palavra código '! 
mais provável, dado a palavra recebida. Serã mostrado que o procedi 


mento descrito a seguir se apresenta mais simples em termos de com- 


. 
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plexidade do decodificador. 


RP E a Considerações Introdutórias 


Em princípio, uma mudança na notação serã realizada por 
conveniência, de modo a torná-la adequada ao problema. Considere R 


o espaço de observações e uma partição de R denotada por to) ' 
je (0,1) e se S, onde S é um conjunto de Índices. É =Iq U n ! 
denota uma classe de confiabilidade, a qual são associados dois pa- 


TÂMB EROd Bo e E 0% Bs. «E Bo < 1; e a quantidade a =Bo- cao 
s s s s s 


GPs s a a ; 
é dita ser o peso da classe Il . Também será considerado um vetor ! 


a! om qu Re ca definido pelo peso das classes nas 
quais as componentes da palavra recebida r = (erre erro) se en- 
s 
contram; ou seja, se re I » então a tl) Rubia Assim, a nota- 
0 


ção (i) indica que se faz referência a classe de confiabilidade a 
qual r; pertence. 

Uma das perguntas importantes que surgem quando é dada 
uma partição (ns) do espaço de observação R, ê como atribuir valo- 


res aos parâmetros Bo e Bo para cada classe de confiabilidade! 
s s 


nº. É claro que o emprego da teoria da decisão seria oportuno aqui, 
usando apropriadamente o conceito de função utilidade [13]. Contudo, 
serã proposta uma maneira de se atribuir valores para os pesos ! 
As de cada classe IÊ da partição, de forma a representarem um 


grau de confiabilidade associado à classe: 


2.5.3.2 - Sobre Atribuição de Pesos a Classes de 


Confiabilidade 


Dada uma partição qualquer (1º), é proposto associar a 


cada elemento Nº da partição os seguintes parâmetros: 
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classe de confiabilidade né => 


g. = pídec correto|r.ellfje 8 =pídec com erro|r elºh. 
é: - e E 


Deve ser observado que 0 <Bo s Bo < 1, e que esta ' 
s s 
atribuição implica em B. + BB, =1. 
ste “a 
De acordo com a definição 2.5.1b, a distância generaliza 


da de Forney entre a palavra recebida r e uma palavra código £ '! 


serã expressa por 


n 
Dç(x,£) = R agr, r£;) (2-31) 
i=l 
onde, 
a s s 
Bo = pídec correto |r,ell k r;7£,+ rel 
Agtrpf)= (2=32) 
o s s 
Be, = pídec com erro|r,ell Jr,7£, rel 


Inicialmente será analisado o caso binário, onde é assu- 
à ESSA i s - 
mido que a partição é simétrica, ou seja, (Il ) é escolhida de tal ! 


maneira que 
pídec O |nº,0) = pídec 1/1º,1) (2-33) 


Para uma classe e, o coeficiente de confiança condicio 
nal de Kiefer |21| é dado por Ry = pídec 3/n82,50. 
Neste caso particular, (2-32) pode ser reescrita na for- 


ma 


ca 0 
Bo = Pídec 3/1" ,5) - rix£,, rjenê 
dgtr;1f))= s (2-34) 


à 0 0 
1-p(dec 3 |nÊ a RA r,enê 


To 
] 


= so " E 
Então para a classe Il |, o peso associado Eq é dado '! 


ea 
por 0.9 = 2.pídec 5] ,9) = 1 =2 Ro T 1 (2-35) 


Notando que 0.5 < Ro gs 1, têm-se que 0 <a “Lo 


Por outro lado, 


eo 
Rjo = ptdec j|l ,jj= (2-36) 


ou seja, 


0 
” pídec 31º |5) 
pidae q Rm eia 
pídec j,11ºº|5) + pídec 1-3,1º2|3) 


Desta forma, deve ser assumido que 


1 E 


Ba = =: — (2-38) 
ºs 1+ Ppídec Vos 188 ado És pit; 4 ]3) 
ptdec 3 ,1ºº]31 ping? [5) 


Definindo Pels" a razão de verossimilhança na classe! 
so 


A png 3) 
Mal O e (2-39) 
elã pts 5) 
e observando que para partição simêtrica fo = fo = 7, a equa - 


ção (2-38) pode ser escrita como 


1 


RE e SR 
Pes Ea 4 consequentemente, ga At TT Tr 


de modo que o peso associado a classe 8º será dado por 
g 

-——— = 2-40 
7 p ( ) 


É interessante observar que o algoritmo de Hartmann - 
Rudolph, descrito no capítulo III, utiliza esta quantidade »p no 
processo de decodificação. 


A generalização deste procedimento para o caso multini 


EsRE 


vel e para casos onde a partição do espaço de observações não é si- 
métrica serã realizada a seguir. Aqui são supostas conhecidas as 


probabilidades a priori dos símbolos da fonte, denotadas por E 


g-1 
testo : 
A equação (2-39) resulta em 
te (a |1º,5) ZE is 
B =.L, Pp. pídec a PS e ae TEM 
RO) ES Es (2-41) 
detr;rf;)= 
1 q-1 é Se si s 
os = “são P; pídec 5/1 ,5) r,=£, rel 
e o peso da classe n8º serã agora determinado de acordo com 
ar Co q 
Gan 2.5£o P; pídec j|i ,j;-1 (2-42) 
g-1 2 O a aca 
+ sho yo P tm 317,5) -5E9P5 (2-43) 
q-1 
= py [2 pídec 5 |n8º,5)- 1] (2-44) 
5=0 


No caso multinível,a equação correspondente a (2-36) é 


O js 
so ptISº |3) 
pídec 5|I" ,5) Sl (2-45) 


q s0,. 
pEoP tp 15) 


= Laden du (2-46) 


qr1 pingo |5) 
ud CEE 


Relembrando da definição de Pp |3 dada pela equação ! 


(2-39), tem-se que 


1 
pídec s(n80,4)= = > (2-47) 
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Desta forma, o termo entre colchetes na equação (2-44)' 


Sitios 3 gs md SE (2-48) 


Da 
e megrá, 
à Sosa a (2-49) 
(e ph 
z é 
?=0 2. |3 


usando o fato que f. = 1, (2-49) pode ser dada por 


5l5 
De Oy 
SR 
Inte sr, = Ts a PO. O ra) (2-50) 
Rr fais 


Esta expressão será definida como Ps multinível, 


ql 
Psi 
a ba as (2-51) 
J g-1 
1 + E ; 
2=0 Fe l5 
t+ 


Então, de acordo com (2-44) o peso da classe nº? será! 


calculado pela fôrmula 


Ela, A P. b (4-52) 


No caso binário, Eq 


de modo que quando a partição ê não simétrica, tem-se Pg * 0 + 


. 
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como exemplificado na figura 2.13 a título de ilustração. No caso '! 


onde (2-33) se verifica, o | = Flo = 9 e consequentemente ! 


Po  Py pi portanto, (2-52) se reduz a (2-40). 


7 


So So 


Figura 2.13 - PARTIÇÃO NÃO SIMÉTRICA, 


E A s F : 

Dada uma partição Il'= (II ), o procedimento aqui estabe 

lecido indica como realizar uma atribuição para os pesos ay das 
classes de confiabilidade, sem indicar contudo como proceder para 


realizar uma "boa" escolha da partição. 


2.5.3.3 - O Algoritmo de Decodificação 


O decodificador proposto deve operar sobre a palavra re 


cebida r, escolhendo a palavra código £ mais próxima segundo o 


critério da distância generalizada, ie, deve escolher uma palavra ! 


código fe C tal que (V£' e C) Dçlr,f) < Dç(x,£") (2-53) 
Para atingir este fato, é proposta uma regra de decodi- 


ficação que consiste em estimar a palavra transmitida como sendo a 


(e) 


= a 
palavra código f e Ca qual minimiza a expressão a(x,E)=Ir (Q) £r.a 


. 


- aja 
Isto pode ser facilmente demonstrado atravês do seguinte lema: 


Lema 1: Dado a palavra recebida r, a palavra código f mais prôxi 
ma segundo o critério de distância generalizada ê aquela que minimi 


za A(r,£f). 


Prova - Dada a palavra recebida r e uma palavra código É, a dis- 


tância genralizada entre elas ê& dada por 


n 
DtEdi SE Ate ot (2-54) 

Se & observado que d.(r.,£f.) = (E. 2E,) ati), gb) 
q Essa dj (x, cf; c ' 


então, & possivel reescrever a equação (2-54) como 


E (1), a(i) 
Dt E) = 4a (rirf,)a Ro (2=55) 
Por outro lado, o segundo termo do somatório independe! 


da escolha da palavra código, de modo que 


n 
- ' (i) 
Min D.(r,£f) => Min E lt à (2-56) 
£ G =—'= £ fez RE 
Finalmente, observando-se que 
Edo (1) 
[r O £].a =; dglrprf,la o, a prova é completada. 
Es02D 


É interessante notar que se existe somente uma classe 


de confiabilidade, então 


Min Dç(x,£) => Min Da (x, É). 
£ £ 
ç 
Pode ser observado que o decodificador primeiramente ! 


efetua uma decisão abrupta determinando uma palavra recebida r, a 


qual pode ou não ser uma palavra código. Em seguida, ele determina 


e 


um vetor de confiabilidade at h, associado a esta palavra recebida. 


É facilmente verificado pelo lema 2 abaixo que se a pa- 
lavra recebida r resultante do uso de decisão abrupta for uma pa- 
lavra código, ela prôpria é a palavra código mais próxima da pala - 
vra recebida no sentido de distância generalizada. Segue-se, portan 
to, que neste caso, o algoritmo resulta em assumir que esta foi a 
palavra transmitida, o que equivale a desprezar as informações de 
confiabilidade. No caso contrário, estas informações são realmente! 
utilizadas no processo de decodificação, possibilitando a correção! 
de erros. 

Lema 2: Se a palavra recebida r é uma palavra código, então a pa- 
lavra código mais próxima de r no sentido de distância generaliza 


da ê a própria palavra recebida. 


Prova - Foi provado anteriormente (Lema 1) que minimizar Do tz, £) é 
o mesmo que minimizar A(r,f), através da escolha de uma palavra cô 


digo, ou seja, 


Min D.(r,£) <=> Min A(r,£) (2-57) 
ad it o 


É claro que 


n 
(vfec) A(r,£) = E ddr f)a?! >0 (2-58) 


<l d= 1,23 csszha 


A conclusão da demonstração decorre do fato que se ea 
& palavra código, então 
Min A(x,£) = Alr,r) = Q (2.59) 
f 


C.o.D 
O lema 2 permite uma simplificação no procedimento de 


=. 


decodificação. Inicialmente uma decisão abrupta é realizada, resul- 


E 


tando em uma palavra recebida r, ea sindrome correspondente a es- 

ta palavra recebida é calculada. 

1 - Se esta síndrome é nula, o decodificador libera a palavra, admi 
tindo que não ocorreu nenhum erro. 


2 - No caso contrário, um vetor de confiabilidade ê& calculado e o 


algorítmo atua da maneira jã discutida. 


2.5.3.4 - Exemplo: 


Serã apresentada uma aplicação do algorítmo de decodifi 
cação descrito, utilizando a treliça associada ao código de bloco " 
binário considerado. Admitindo que a escolha dos pesos das classes! 
de confiabilidade utilizadas esteja de acordo com o que foi propos- 
to (2.5.3.2), e que a partição simêtrica foi escolhida de modo a 


resultar em uma máxima variabilidade dos pesos das classes; a regra 


de decodificação pode ser anunciada como: "Estime a palavra transmi 


n : 
tida como a palavra código para a qual 181 dtg, ef, o 8) é mínimo, 
onde, p (É) =1- qi) 1 + gti) e g 6) é a razão de verossimi - 
lhança”. 


Este procedimento pode ser facilmente implementado com 
o auxílio de uma treliça. No exemplo, novamente é assumido que (o) 
código (5,3,2) é empregado, e que o canal é perturbado por um ruído 
com distribuição de probabilidade N(0,1). É admitido ainda que a ! 


fonte bipolar transmite +1V ou -1V correspondendo aos símbolos ! 


0 el. 
Suponha que as amostras recebidas foram ' 
(-0.7,-0.2,-0.3,+0.7,+1.1),de modo que pf s- (0:76 ,.052,5043,0,6,/0,8] 
A palavra recebida ê aquela que seria decodificada por 
decisão abrupta, ou seja, r = (0,0,0,1,1). 


Usando a treliça, são obtidos os seguintes passos no 


processo de decodificação: 


E 


0.0 0.0 0.0 
0.8 0.3 

0.2 

0.6 0.5 


(a) Vad tás 


0.2 
0.5 


(d) 


Figura 2.14 - PASSOS NA DECODIFICAÇÃO POR DISTÂNCIA 
MÍNIMA GENERALIZADA. 


A palavra código (1,0,0,1,1) estã mais próxima da pala - 
vra recebida (0,0,0,1,1) no sentido de distância de Hamming; entre- 
tanto, a palavra decodificada foi (0,1,1,1,1), mais prôxima no senti 
do de distância generalizada. Os dígitos corrigidos foram exatamen- 
te as posições de mais baixa confiabilidade. Tambêm deve ser obser- 
vado que a escolha da palavra código (1,0,0,1,1) implicaria em in - 
verter um dígito com confiabilidade relativamente alta. Finalmente, 
este procedimento apresenta uma vantagem computacional, pois nem 
sempre é necessário efetuar cálculos, visto que quando dA (x 0f,) =0, 


nada se adiciona ao somatório. 


2.6 - DECODIFICAÇÃO POR DECISÃO SUAVE EMPREGANDO SUBTRELIÇAS 


Matis e Modestino |25|] introduziram têcnicas que permitem 
uma redução na complexidade da decodificação de códigos de bloco li- 
1 


neares para algoritmos que empregam a treliça gerada pela matriz 


[H]. Uma das formas proposta para realizar isto ê atravês do uso de 
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subtreliças bem mais simples que a treliça associada ao código. Se 


a palavra recebida r tem vários coeficientes de confiabilidade ! 


com valores grandes, as saídas de decisão abrupta correspondentes '! 
podem ser consideradas corretas. A subtreliça é gerada aceitando co 
mo corretas estas saídas de decisão abrupta para os símbolos nas '! 
k-p posições mais confiáveis, e permitindo todas as possíveis com- 
binações nas posições restantes. É interessante observar que a esco 
lha p=k corresponde ao emprego total da treliça, enquanto que (o) 
extremo oposto p=0 corresponde a realizar uma decisão abrupta sem 
uso da treliça. Valores intermediários de p fornecem as subtreli- 
ças desejadas. A aplicação deste procedimento é melhor ilustrada ! 
considerando o exemplo apresentado na sec. 2.5.3.4, onde 


r = 10,0,0:1,11 cul = (66,02,03,06,08). 


As subtreliças geradas para este caso são mostradas na 


figura 2.15 (a) e (c). 


(a) P=2 


(b) Pp=1 (co) p=1 


Figura 2.15 - EXEMPLO DE SUBTRELIÇAS PARA UM 
CÓDIGO DE BLOCO. 
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O desempenho da decodificação por máxima verossimilhan- 
ça (descrito na sec. 2.3) empregando subtreliças ê analisado para 


alguns côdigos, na referência |25|, através de simulação em computa- 


dor. 


CAPÍTULO III 


DECODIFICAÇÃO ÓTIMA DE CÓDIGOS LINEARES 


Serão apresentadas duas técnicas ôtimas de decodificação 
de códigos lineares em canais discretos no tempo e sem memória, quan 
do as palavras código são equiprováveis. Ambas fazem uso de decisão! 
suave e são ôtimas no sentido que minimizam a probabilidade de erro 
por símbolo e por bloco, respectivamente. A primeira delas, o algo - 
ritmo de Hartmann e Rudolph [18], faz uso das propriedades de linea- 
ridade do código, sendo portanto, aplicável somente para códigos li- 
neares. Ela é particularmente atrativa para côdigos corretores de '! 
erros com alta eficiência, conforme serã visto. A segunda técnica '! 
descrita é um algorítmo para maximização da probabilidade a posterio 
ri das palavras código |29], aplicável a côdigos não lineares tão 
bem quanto para códigos lineares. Quando estes são considerados, fo) 
uso da treliça associada na decodificação torna prático seu emprego! 
por simplificar o processo de decodificação, preservando ao mesmo 
tempo sua otimalidade. O desempenho deste algoritmo & inferior ao 
desempenho da regra de decodificação de Hartmann-Rudolph com relação 
a taxa de erros por símbolo, e vice-versa com relação a taxa de er - 
ros por palavra. Para ambos os algoritmos serã analisado o desempenho 
assintótico da probabilidade de erro, por bit e por palavra, respec- 
tivamente, quando são utilizados códigos de bloco binários e o canal 


é perturbado por ruído branco gaussiano. 
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3.1 - O ALGORÍTMO DE HARTMANN E RUDOLPH 


3.1.1 - Introdução 


A regra de decodificação apresentada pelo algoritmo de 
Hartmann-Rudolph ê ôtima, no sentido de que minimiza a probabilidade 
de erro por símbolo sobre um canal discreto no tempo e sem memória , 
para qualquer código linear corretor de erros ,quando as palavras cô- 
digo são equiprováveis. 

A maioria das têcnicas de decodificação de códigos correto- 
res de erros são exaustivas no sentido de que toda palavra código & 
utilizada no processo de decodificação. Estas técnicas não fazem 
qualquer uso essencial da linearidade do côdigo: O algoritmo aqui 
apresentado ê também exaustivo, mas no sentido de que toda palavra 
do código dual é utilizada no processo de decodificação; por isto |, 
na prática, esta regra ê usada somente para códigos cujo côdigo dual 
tem um pequeno número de palavras, ie, ela ê particularmente atrati- 
va para códigos de alta eficiência. Por exemplo, para o código de 
Hamming binário C(15,11,3) existem ah palavras código, enquanto 
que o código dual C'(15,4,3) possui apenas as palavras código. É 
concluido, portanto, que a complexidade da regra de decodificação de 
verá variar inversamente com a eficiência do código. Este algoritmo! 
é de natureza essencialmente estatística e, é aplicável tanto a códi 
gos de bloco quanto a códigos convolucionais. No apêndice A ê apre - 
sentado um programa em Fortran IV utilizado na investigação do desem 
penho deste algorítmo para códigos de bloco binários em canais com 
ruído branco aditivo gaussiano. Os resultados destas simulações se - 


rão discutidos no capítulo a seguir. 


* Uma maneira de utilizar a linearidade é considerar o uso da treli- 


ça associada ao código, descrita na sec.2.3. 


mOs 


3.1.2 - A Regra de Decodificação 


A regra de decodificação ê apresentada para códigos de 
bloco lineares, e posteriormente uma extensão para códigos convolu- 
cionais serã feita. 

Seja c= (cgrSyrece1Spo9) uma palavra de um código de 


: t = , , 1 
bloco linear C(n,k,d) sobre GrF(p), e £j Tas Rodada oe OUR 
j-êsima palavra do código dual c'(n,n-k,d'). A palavra c é J 
transmitida atravês de um canal discreto no tempo e sem memória cu- 
jos símbolos Es são números reais, de modo que a palavra recebida 


é da forma r = (rg'ryr-*e+En0))+ com r, e R. A questão que o algo 


ritmo se propõe a resolver ê: dado r, calcular uma estimativa Em E 
do m-êsimo dígito Cm da palavra transmitida c, de maneira que 
a probabilidade de se ter E igual a ca seja maximizada. De ou- 


tra maneira, a estimativa êm é tal que 
pule = Em lx) > Pr(c) = sr), vs e GF(p). (3-1) 


A regra de decodificação pode ser assim descrita: Esti- 
me o dígito transmitido c, como sendo o valor de s e GF(p) que 


maximiza a expressão 


n-k 


p= 2 n-1 p-1 -ilc!,- ) 
Agts)=2 0 ace TSE CO jet Pr (rp li) |; (3-2) 
t=0 j=1 [£=0 i=0 
ou seja, a estimativa e deve ser tal que An (Em) = Max Ants). 
seGF (p) 


Isto pode ser provado mostrando que P(co=s |) =A Ants) 
onde À & uma constante positiva. 


Inicialmente, ê possível escrever que 


pr(c =sl|r) = LL Plclr) (3-3) 
= L CE. 
ec, 
Css SE 


- 


AQE 


ou 
pr(lc=s|r) =.z P(rlc) [P (2) /P(1)] (3-4) 
ea ES ca E 
cel, 
c=s 
m 
Considerando que os símbolos de C são elementos de 


cr(p) e que as palavras código são equiprováveis, tem-se que 


P(c) = see, e (3-4) pode ser reescrita como 


(3=5) 


Onde 825 é o delta de Kronecker, 


a se L=j5 
Ns 


e e 4 (8 mo 
0, em caso contrário 


“..1Ô 
=m É sm 


ml n-1 


Em termos de transformada finita de Fourier, pode ser ! 


mostrado que |18| 


pel 
tam” RP E coa (3-6) 
t=0 
e 
o a 
P(rlco =P” 2 F(r,u)6 (3=7) 
uevV 
ba 
=U.V 
onde F(r,u) = E P(rlv)o (3-8) 
“Va 


e 0=exp(v-1 27/p) representa a p-ésima raiz complexa da unida- 
de, e V5 O espaço n-dimensional sobre GF(P). 
Substituindo (3-6) e (3-7) em (3-5) e levando em conta 


que devido as propriedades de ortogonalidade de caracteres de grupo 


rd E Be mise! 
ÉS = 
1» (3-9) 


0 , em caso contrário 


“ole 


É obtida a expressão 


-n-1]p-1 e ares 
Pie selr) = |E-—| 2 6 LT FPl(r,c!- te) (3-10) 
Rasa P(X) Jt=o 31. Cd cm 


Para canais sem memória, (3-8) pode ser reescrita como 


agi -u, v, Rol pel -iu 
F(x,u) = E n Plr,|vp)o tem E P(rp |i)0 t ê de 
vev t=0 £=0 i=0 
modo que se obtêm a expressão 
-n-1 PH é 
P(c =slr) a RREO 
m Pl!  E=0 
oe el usigdo otpgetá 
es RM a o) EVA me p(rpli) (3-11) 
j= = ts 


Definindo A = pr Lt), e observando a definição de 

Ants), segue-se que P(c =slr) = À Ants) 
C.o.D. 
É interessante fazer uma aplicação deste algorítmo no 
caso de códigos lineares binários, onde a regra de decodificação se 
apresenta de uma forma comparativamente mais simples. Neste caso, 
p = 2, e o procedimento para decodificação consiste em escolher a 


estimativa 


0, se AntO) > An (1) 
fo) - (3-12) 


ls em caso contrário 


Esta forma comparativa da regra de decisão pode ser a 
presentada de maneira mais conveniente, quando estabelecida em ter- 
mos da razão de verossimilhança 9 = Pr (ro |1)/Pr(r 0). 


. 
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Para o caso em questão, a desigualdade Ano) > An (1) re 


sulta Ee 
12 -d(ol, = +84] 
Ro NT Do qr J no pr(rpli) > 
t=0 j=1 2=0 1=0 
1 ge 4 -i(cl, = t5 
E EINS I E (=1) j * petrç]t) (3-13) 
t=0 j=1 2=0 1=0 
ou seja, 
22"kn-1 cl, - 8 
z E ircte,/0) + tem * ne prtr,[1))> 0 43-14) 
j=1 2=0 


Dividindo ambos os membros de (3-14) pela quantidade po- 
nel 
sitiva I Pr(rp|0) e usando a definição de razão de verossimilhan 


ça, a desigualdade pode ser escrita como 


29"k nu2 RR 
z m Ê +9, (1) O nt) > 0 (3-15) 
j=1 450 


Agora dividindo (3-15) pela quantidade positiva ú 


n=l 
TP ql Pp) e utilizando a identidade 
L=0 
-Cip - Sm cie O mt 
1 + pato = 1-9, 
= |[-—— + Segue-se que 
1 + Ap 1+ %p 


a regra de decodificação no caso binário pode ser estabelecida con- 


forme enunciada abaixo: 


E cl,(0)6 
9 k SE, aa % JL mt 
Escolha & =0Q se 1 1 |——— >» O; “Caso 
A JEM p=oMaa 
PA 
contrário faça ce aum 


Algumas observações devem ser feitas. Em primeiro lugar" 
deve ser notado que quando esta regra de decodificação simbolo-por- 


simbolo ê usada, a palavra estimada & na saída do decodificador ! 


A 
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não é necessariamente uma palavra código. Em segundo lugar, se côdi 
gos de bloco cíclicos são considerados,a regra de decodificação de- 
ve ser determinada para o simbolo recebido r9r € então, os simbo - 


los restantes Free l og podem ser decodificados simplesmente '! 


permutando ciclicamente a palavra r no registro a deslocamento '! 


(acumulador), como ilustrado no exemplo a seguir. 


3.1.3 - Exemplos 


Serã utilizado o código de Hamming (7,4,3) para exempli 
ficar a regra de decodificação binária aplicada a códigos de bloco, 


que neste caso torna-se: 


cir O Sme 


8 6 1 - Mp 
Escolha q =0Q se £ 1 |-——— SO), e 
j=1 2=0 | + Pp 


e q = 1 em caso contrário. 


A matriz de verificação de paridade deste código é 


À o o Cais A 
H]=|0 1 da E 2d 
PS eo; vd. + 


e consequentemente o espaço linha c' serã dado por 


sido 


Cc! 


0 T 0 0 1 1 1 
1 0 a 0 0 1 J 


Figura 3.1 - ESPAÇO LINHA DO CÓDIGO DUAL DO CÓDIGO DE 
HAMMING (7,4,3). 


Tomando Pp ê (1 - Fp/a + Pp) a regra de decodifica- 


ção consiste em escolher & = 0 somente se 


P9t'PyPaPa + PoPsl6 + age Lat E PaP4Ps + PoPrP2º3º5 E PoPaP3º4P6* 


4 BRAD atada * Os 


O decodificador correspondente estã esquematizado na 


figura 3.2 exibida abaixo 


oe 1 


1 caso contrário 


Figura 3.2 - DECODIFICADOR ÓTIMO PARA O CÓDIGO DE BLOCO(7,4,3). 


. 
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Para ilustrar a aplicação desta regra de decodificação! 
para códigos convolucionais será usado um código convolucional E 
(4,3,3). A regra de decodificação para o símbolo recebido r9 usan 
do um código convolucional (ng rkgrN) será, no caso binário, enuncia 
da como: 
' 
o 4a e; ciç O me 


Escolha q = 0 se e somente se £ 1 |——— > 
j=1 t=0l1 + 9, 


de outra forma assuma So =-Ls 


Naturalmente, existe somente um número finito de termos 
não nulos na equação acima, dependendo do comprimento da sequência" 
transmitida. 

Agora, para exemplificar, as porções iniciais da matriz 
de verificação de paridade descrita no exemplo 1.5 e do espaço li - 


nha Cc! são mostradas abaixo. 


+ 
O Msosmntsadanda ca das po o tom dn dns d 0 sa 66 6/0166 “ 
ART cho DE si Aa aa ceras e alga 666 Seca VE 65 
é E DD TS O secos ese sa 


QI O DLL Oparsanty ais qe masa 


Dil, COMO ado Ci o ROLO O Iso ietçia]s 
+ + 
Figura 3.3 - ESPAÇO LINHA DO CÓDIGO DUAL DO CÓDIGO 
CONVOLUCIONAL (4,3,3). 
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Como antes, pp A A - Ap1/( + Ap) de modo que a re -— 


gra de decodificação serã expressa por: 


ê&y = O se e somente se py + P1P283 + PoPaPsP6Ps É PoPyPaPalsPofg + 
t «e. > O, em caso contrário faça Eç-1. 

O diagrama do decodificador & então mostrado na figura" 

3.4, tomando forma de diagrama de treliça para o código dual E 
c'(4,1,3) com as posições Cio nos rôtulos das malhas complementa- 
das”. (Em geral, para decodificar rp! às posições Cim devem ser 


a, 


complementadas). 


PotPIP2P3 (P5+P9P 123) PaPsP603 +lpgteyPaP3)1 


SoPifatea 
PorPyrPorhs Pato agud 


Figura 3.4 - DECODIFICADOR ÓTIMO PARA O CÓDIGO CONVOLUCIONAL 
(AS 3) 

Desde que diferentes unidades de memória devem ser em - 
pregadas para cada símbolo a ser decodificado, a quantidade de memô 
ria requerida por este tipo de decodificador cresce linearmente com 
o comprimento da sequência transmitida. O mesmo também é verdade pa 
' 


ra o algoritmo de Viterbi aplicado para decodificação de códigos 


convolucionais. 


* Veja figura 1.4, diagrama de treliça para o código convolucional'! 


(4,1,3). 


SAT 
3.1.4 - Desempenho Assintótico do Algoritmo 
Expressões assintôticas para a probabilidade de erro por 


bit na decodificação de códigos de bloco binários lineares em canal" 


com ruído branco aditivo gaussiano fóram demonstradas |19 


A regra de decodificação ótima bit-a-bit para o m-êsimo'" 
dígito de um código C(n,k,d) linear e binário consiste em escolher 
Em = s, onde se GF(2) e maximiza PE, = s|r). Desta forma a pro- 
babilidade de erro por bit é dada por 


Ta P[P (Em = Cmlz) g P(Er Cd Six) | (3-16) 

Não hã perda de generalidade em assumir que a palavra 
transmitida foi a ênupla toda nula, desde que o código É linear e o 
ruído apresenta simetria. Para facilitar a dedução, é considerado 
que a palavra toda nula foi transmitida. Quando c = 0 é transmiti- 
da, a m-ésima componente da palavra recebida r é E VE + E ai 
onde E é a energia do sinal por bit e em uma amostra de ruído de 
um processo gaussiano com densidade espectral de potência unilateral 
no watts/hertz. A SNR para este canal é Yy = E/nogr ou em termos dos 
bits de informação transmitidos, ia Es/no RÃ EP e m-êsima 


componente da palavra recebida r será decodificada incorretamente! 


se e somente se 


plc, = 0lz) < Plc = 1]x) (3-17) 


onde fis (VE + egre cc rYE + e 


agilia 
Em outras palavras, 


Pora = Plato  PtejEilas =p. Etelt) (3-18) 
bt (eis, E AR CES) | 


. 
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onde É, A fo e Cleo = 4), = Oia 
a = m 


Desde que o canal é sem memória e as palavras código são 


equiprováveis, (3-18) pode ser escrita como 
Po ra RU I P(rplcp) E I P(rplcp) (3-19) 


c 
Como rp =ep + (=1) t VE & uma variável aleatória real 


e contínua com densidade de probabilidade 


o 
(xp = (10) 4/5) 


END RCE, TED e Spa 1 
g(rplcp) exp 7 e o produto 
Mo 0 
n-1 
IH g(rp|0) é positivo, segue-se que 
L=0 
n=1 g(r,|€,) n-1 g(rplc,) 
PiçieBls pn Eder q e (3-20) 
cesy t=0 g(r,|0) ces, £=0 g(r,10) 
ou seja, 
Pasto P] É exp (-r.cvRY,/n9) < z exp (-r.cvRYy/9) (3-21) 
ces ces 
Aloe 1 
Para análise do comportamento assintótico de Pot” deno 
e “ 1 1 
tado AB, p+0 Phit! e AB tre Pripl é no será fixado e E, variado 


de modo que VE varie. 


3.1.4.1 - Caso de Baixa Relação: Sinal-Ruído 


Para x << 1, exp x * 1 + x, de modo que numa pequena vi 


zinhança de zero é possível escrever 
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AB, 50 Ppit)z “e: (1 = r.c 4/RYp/NQ) < 


ces) 
7 O-re 4 'o| (3-22) 
ces 
="1 
Então, como |So| = |s1 5 


AB, p»0(Ppit) =P E = .E.e 4/RYp/No< 20 ap eis 4/RY p/o | (3=23) 


Desde que s5 & um código linear binário (n,k-1), se - 
gue que se os vetores de So são arranjados como linhas de uma ma- 


= ” k-2 
triz Mg! então cada coluna serã toda nula ou conterã 2 zeros 


e gra l's. Arranjando o conjunto de vetores de s, como linhas 
de uma matriz Mr as colunas de M que correspondam a colunas to 
da nula em Mg serão toda um, enquanto que as demais colunas terão 
io zeros e ias EB 
Usando este fato, é possível escrever 
Jo 
BB po (Ppit) = "fe; rp £ | (3-24) 
31 


onde dprdoreccrdo são colunas nulas de Mp - 
Desde que rp - N0E:nç/2) para L=3y.ccgr 
+oo 


P| EL r,;<0|=—=—— exp(-x2/2) dx (3-25) 


0RYj 


De (3-24) e (3-25) segue-se a expressão assintótica de- 
sejada, AB, p»0 pit) = Q(/Z6RY) - 
Se o código dual de C tem distância mínima maior que 


2, então 0 = 1, e o comportamento assintótico pode ser descrito pe 
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la equação AB, pq (Ppit)z Q('2RYp). 


3.1.4.2 - Caso de Alta Relação Sinal-Ruído 


Desde que 0 € So! Z exp(-r.ce 4/RY6/No) > à RR 
ces 
=""0 


Se r.c > 0 para todo "cEs,, então quando + o 


si + 


PjI<i exp(-r.c 4/RY,/No) <L expl-r.e 4/RYp/n0) > 0 (3-26) 
cesg ces 


a 


consequentemente, (3-21) é cotada superiormente pelo "union bound " 


quando Yp Cresce, ie, 
AB (Prsal sobf race seo] A oie 20) (3-=27) 
ce É 


Pode ser mostrado |19| que esta cota é apertada para va 


lores de Y suficientemente grandes, de modo que 


b 
AB Polo Ro Bixo <0) (3-28) 
Yyb-+ bit es si ei om 
E 
Assim, +oo 
Pes oue E exp(-x2/2) 6x = Q(V2RW (Cc) Y,) 

Soda y21 (3-29) 

Y2RW(C) Yy 


Portanto, substituindo (3-29) em (3-28), tem-se; 


: p= -—- 

AB pro Ppit) = E  Q(W2R W(c)Yy) (3-30) 
CES 
me SE É 

Se Wm é o menor peso das palavras código em sj e 


N(W) o número de palavras côdigo de peso Wm! então 


Odo. 


( )= N(W) Q(/2R Wo Yo) isto porque para valores suficien- 


Phit 


temente grandes de x, Q(x)z 1/xv2I exp(-x2/2), de maneira que so - 


AB oo 


mente as palavras côdigo de peso no contribuem de forma signifi - 
cante no somatório. Deve ser também observado que se o côdigo é ci- 


clico e tem distância mínima d, então “E =d e 


BB, jo o Ppip)= N(d) Q(V2Rd Yy) 


3.2 - ALGORÍTMO PARA MAXIMIZAÇÃO DA PROBABILIDADE A POSTERIORI 


3.2.1 - Introdução 


Uma outra regra de decodificação serã a seguir apresen- 
tada a qual é tambêm Otima, mas no sentido de que para qualquer cô- 
digo corretor de erros minimiza a probabilidade de erro por palavra 
sobre qualquer canal discreto no tempo e sem memória, quando as pa- 
lavras código são equiprováveis. Esta regra implica em maximizar a 
probabilidade a posteriori das palavras côdigo. É sabido que exceto 
para poucos canais clássicos não é normalmente muito claro como de- 
ve ser processada a palavra recebida de modo a maximizar a probabi- 
lidade a posteriori. No que segue é estabelecido um procedimento pa 
ra realizar tal maximização. Deve ser observado que ao contrário do 
algoritmo de Hartmann-Rudolph, a decodificação resulta sempre em 
uma palavra código quando este procedimento é empregado. Serã tam - 
bém mostrado que a regra de decodificação para o caso binário pode 
ser enunciada como: 

Escolher a palavra código c que maximiza a expressão! 

E cptn Pp 


Se dado a palavra recebida r = (rgrr ) for 


qe 


* 
definido um vetor r calculado a partir de r; 


* 
r = (Log Agrtog P,r--«ntog Pra onde 9 ,=Pr (x, |1)/Pr (x. |0), o 


. 
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decodificador deve escolher a palavra código -c que maximiza A 


* * * 
c.r , ie, &r = Max C.r , Isto torna possível uma interpretação in 
(e) 


teressante, observando que este procedimento generaliza a decodifi- 
cação por correlação para códigos de bloco sobre canais com ruído ! 
aditivo gaussiano, uma das técnicas Ôtimas conhecidas no sentido '! 
de que minimiza a probabilidade de erro por bloco quando as pala - 
vras código são equiprováveis. Deste modo, é plausível interpretar"! 
esta decodificação como sendo uma forma de "correlação generalizada 

Quando códigos lineares são usados, o emprego da decodi 
ficação usando treliça faz este receptor interessante por simplifi- 
car as operações de decodificação, retendo entretanto, a otimalida- 
de. É importante observar que este algoritmo de decodificação é 
atrativo tanto para códigos com alta eficiência, quanto para os de 
baixa eficiência. Se o código utilizado ê de baixa eficiência, en - 
tão a decodificação pode ser feita exaustivamente, ie. utilizando 
todas palavras código no processo de decodificação. Ele também e 
de particular utilidade na decodificação de códigos com altas taxas 
visto que a complexidade da treliça é função do número de dígitos ! 
de paridade. 

Se ao invês de empregado como na forma descrita acima ! 
o algorítmo for utilizado quantizando em 2º regiões as amostras '! 
recebidas, o resultado é uma pequena degradação no desempenho, po - 
rêm, facilitando sobremaneira a implementação do decodificador, co- 
mo mostrado a seguir. A decodificação agora ê feita atravês do se - 


guinte procedimento: 


1 - vetor recebido Eae GegrryreceXo) com r, E R 
- 1 = 1 , , 1 
2 - vetor quantizado r (rgrrjreccerçoo) com rj € Ro 
3 - vetor Log razão de verossimilhança 
* * * * * 
x =(xgelyrec er) com r; € Rg 
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R = (espaço das possíveis saídas do canal ruidoso) 
Ro- (espaço das possíveis regiões de quantização) 
Rg- (espaço dos valores do fog razão de verossimilhança) 
É então utilizada a treliça associada ao código para ma 
é *& nel * 
ximizar a "correlação" c.r =L C,T. 
Es 2=0 


O algorítmo serã deduzido para códigos de bloco com al- 
fabeto de simbolos p-ários, supondo palavras código equiprováveis ! 
transmitidas em um canal discreto no tempo e sem memória, e poste - 


riormente uma extensão para códigos convolucionais será feita. 


3.2.2 - A Regra de Decodificação 


Seja c= (CgrSy rec e+S40)) uma palavra côdigo de um có 
digo de bloco C(n,k,d) com símbolos em GF(p). Ao receber a pala - 
Vea = (rgrryrec erro)» o receptor deve decidir pela palavra cô- 
digo c que maximiza a expressão 

n-l p-l pal 5(cp-i) 


mo 2 Plrlid E 
2=0 1=0 j=0 


(3-31) 


onde O = exp|2Il/-1/p| representa a p-êsima raiz complexa da unidade. 
Serã provado que a palavra código que maximiza (3-31) '! 
também maximiza a probabilidade a posteriori P(c|r). A probabilida 


de P(clr) pode ser determinada por 


P(c|lr) = P(c) P(r|c)/P(x) (3-32) 
e como por hipótese as palavras são equiprováveis, P(c) = pe de 
maneira que 
+ ao 
P(c|jr) =P — p(rl|c)/P(x) (3-33) 
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Em termos de transformada finita de Fourier, a probabi- 


lidade P(rlc) pode ser escrita como 


Ptrle) =p” 1 (ruoê'£ (3-34) 
uev 
— nQ 
onde F(r,u) = E P(r|vio EL (3-35) 
VEN a 


Como são considerados apenas canais sem memória, 


n-1 -UpeV 
F(r,u) = E I P(rp |vp)O Ms 
eli vev, t=0 
n-1 prl -iup 
= 1 z P(rp |i)o (3-36) 
t=0 1=0 


Substituindo as expressões (3-36) e (3-34) em (3-33) re 


sulta 


esa n-1 p-l 4plep-i) 
P(clr) = [p /P(0)] E T E P(rp]i)o 
uv, t=0 i=0 


-1 p-1 =1 
psd Po s(cp-i) 

I z P(r, li) E B (3=37) 
=0 1=0 j=0 


1] 
112 


De modo que maximizar (3-31) pela escolha de uma pala - 
vra código & o mesmo que maximizar a probabilidade a posteriori das 


palavras código. 
C.Q.Da. 


Serã feita a seguir uma aplicação deste algoritmo para 
o caso binário, ie, quando p=2. 


n-1 cp cp 
P(clr) = [Pa] Ir |P(rpl0)(1+0 ) + P(xp]1) (1-0) (3-38) 
2 


miQiG 


Mais uma vez é conveniente exprimir os resultados em ! 


termos da razão verossimilhança 7, = Pr (r, |1)/pr (rp |0) 2 e de 
Pp A (1 - Pp) /( + Pp). Isto resulta em 
SA n=1 c c 
P(clr) = [2 x Pyp(r) I E + 0 + + Pp (1 a B.zplo 
ae 
SEE n-1 
Tomando A = [2 /P(x) | P(r,|0), tem-se que 
TO R=0 
n-1 cp 
P(clr) = À : LO + 8) + 0 SA - 7,)] 
n-1 cp 
= À elo 1 + Pp (o) e(l + Pp) (3-39) 
+ se c, =0 
c SET: : id 
Contudo, O t = exp Ae = 
1, se Cj = 1 
fg E 
ou seja, 0 = (-1) + de forma que (3-39) pode ser reescrita como 
n-1 cp 
P(alzi = A To Il +ppaltiã) (1+9,) (3-40) 
e 2=0 
c 
cp 209" A 
Usando a identidade 1 + Pp(-1) =>———— |, obtêm-se'! 
1+9, 
n-1 cp 
P(clr) =2A IN fp , e A independe de c. 
t=0 


Como o logarítmo é uma função monotônica crescente, a 


maximização de P(c|r) no caso binário pode ser realizada escolhen- 


do a palavra código c = (core *-"Ch=9) que maximiza a expressão '! 


n=-1 
Z O End 
De ça 
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3.2.3 - Aplicações para Côdiyos Lineares 


É fácil verificar que no caso da transmissão de sinais! 


polares em um canal com ruído aditivo gaussiano, a regra de decodi- 


: n-1 
ficação consiste em maximizar a quantidade c.r =1L Cp Lp pela 
O  R=0 
escolha de uma palavra código e= (CgrSyre cc Spo))- Se os sinais ! 
transmitidos são pulsos liga-desliga de amplitude /E volts ou 0 


volts, então a cplrp - VE/2) deve ser maximizada pela escolha de 
Ea k 

Numa aplicação particular serã obtida a regra ôtima de 
decodificação quando são empregados detetores de envoltória. Consi- 
dere que o sinal ê transmitido atravês de pulsos de RF liga-desliga 
(ASK) em um canal com ruído aditivo branco gaussiano. É assumido '! 
que a envoltôria da forma de onda recebida & amostrada em interva - 
los bastante distantes de modo que seja razoável supor que as amos- 
tras são estatisticamente independentes. O sinal recebido é 


Z(t) = A, cosw + xa (t) Cosug+ — Xo (t) sen WU, (3-41) 


ot 
onde Ag =(0 e Aj = constante, e o ruído estã escrito na forma de 
banda estreita |33| de modo que xçtt) e x (t) são variáveis alea- 
tórias iid gaussianas. 

A saída r(t) do detetor de envoltória ê dada por 


djs DE LO, (3-42) 


2 2 
= + 
r(t) (a x, 664] + xo (6) 
A densidade de probabilidade da envoltória em um instan 
te de tempo t [33] tem a distribuição de Rayleigh na ausência de 
sinal, e uma distribuição de Rice quando o sinal está presente, ou 


seja, 
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SA 
f 2 2 
= exp (-r,/20 er rp 2 0 
P(rp|0) = (3-43) 
0 , rp <0 
e 
2 2 
ê r, + A Aos 
L Tê sb Eos 
-— exp - I FED 
q2 202 0 52 L 
P(rp|1) = 
0 E rp < 0 (3-44) 


onde a* = Etnô(t)) = EtnÊ (6)) é a potência do ruído, e Ig (x) ê 


a função de Bessel modificada de primeira espécie e de ordem zero . 


P(r 


aba 


Figura 3.5 - DENSIDADES DE PROBABILIDADE DA ENVOLTÓRIA. 


A razão de verossimilhança ê dada pela expressão 


P(rp|1) Dna 2 
Yp = Plrp TO = exp(-A//20 ) Ig (xçAy/0 ) (3-45) 


Definindo a relação sinal/ruído como Yy = n$/202, o va 


lor do limiar b pode ser encontrado resolvendo a equação ! 


exp (-Y) Ig (bg v/27) = 1, cuja solução pode ser obtida atravês da exce- 


lente aproximação |34| 


=) Bus 


by = v2+Y/2=b/0 (3-46) 


Uma análise do comportamento assintôtico da regra de de 
codificação para baixa e alta relação sinal/ruído serã realizada a 


seguir. 


3.2.3.1 - Baixa Relação Sinal/Ruído: y << 1 


Neste caso ê quase sempre verdade que |29| 


v RO 2 
Ig (r,A,/0 j= explry Aj/40 ) (3-47) 
portanto, 
Mo | m 
Rn Xp a is IL (3-48) 
[o] 20 


e a regra de decodificação consiste em escolher a palavra código '! 


c'= (cgrre-1Cn=)) que maximiza a expressão 
n-1 E 
Raio sms (3-49) 
2=0 t 20 


Para ilustrar a aplicação desta regra de decodificação, 
serã considerado o código de bloco binário (3,2,2) cuja treliça as- 


sociada estã apresentada abaixo 


10 


Figura 3.6 - TRELIÇA ASSOCIADA AO CÓDIGO LINEAR (3,2,2). 


A 
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É assumido que a ênupla recebida & r=(1.7,.1,.5), e 
que o ruído tem média zero e variância unitária. Os passos no pro- 
cesso de decodificação estão mostrados na figura 3.7, resultando na 


palavra código (0,0,0). 


ao o | AS TESS, RE! 
445 445 


Figura 3.7 - PASSOS NA DECODIFICAÇÃO PARA O CÓDIGO (3,2,2). 


Deve ser observado que o emprego de decisão abrupta re- 


sultaria na palavra (1,0,0), uma vez que os símbolos rp são compa 


rados com o limiar b = 0/2 + Y/2 = oy2. O erro corrigido foi exata 


(r, = 1.7), ou seja, a 


mente na amostra mais próxima ao limiar 0 


amostra de menor confiabilidade. 


3.2.3.2 - Alta Relação Sinal/Ruído 


Neste caso a distribuição de Rice pode ser aproximada ' 


por uma distribuição gaussiana |6 | 


2 
E A 
1 pé 1 
P(rp|1) == exp - 5) 
2Ilo 


(3-50) 


De modo que o tog da razão de verossimilhança pode ser 


dado aproximadamente por 
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Mm. A». o2 
Enf, É ——pem (rp- ES sé ds E. (tnrp-tnoç) (3-51) 
onde cy = o/v2I. 

Portanto, a palavra código escolhida deve maximizar a 


expressão (3-52) abaixo 


n-1 o 
E c,| (rp=h/2)0 — =(Lnt; = Eno) (3-52) 
28 | Vaga? x 2 o| 


A extensão para aplicação desta regra na decodificação! 
de códigos convolucionais (ng kg) consiste, no caso PinRTaA, em 
escolher o caminho da treliça que maximiza a expressão Eae caia , 
onde novamente o número de termos não nulos depende do comprimento 
da sequência transmitida. 

O algoritmo de Viterbi para decodificação de côdigos 
convolucionais pode ser facilmente modificado para fazer uso desta 
regra de decodificação, e a melhoria provida pelo uso desta técnica 


de decisão suave encontra-se analisada na referência |1 |. 


3.2.4 - Análise do Comportamento Assintótico 


Serã deduzida uma expressão assintótica para a probabi- 
lidade de erro por palavra na decodificação Ótima de códigos de blo 
co lineares binários, usando decisão suave , para canal com ruído ! 
branco aditivo gaussiano, quando a relação sinal/ruído é elevada. A 
expressão obtida é função da distribuição de pesos das palavras cô- 
digo e da relação sinal/ruido. 

Aqui a atenção serã restrita aos códigos de bloco linea 
1 


res binários com palavras côdigo equiprováveis transmitidas sobre 


um canal com ruído branco aditivo gaussiano atravês de pulsos liga- 


-1Dl- 
desliga. É assumido que o código C(n,k,d) considerado tem palavras 


EIA K k 
código denotadas por c., i=0,1,...,2-1, onde Si (Cigregreees 


à 


Cc, j 
in-1 
A regra Ótima de decodificação é escolher & = sec 


1 


que maximiza P(C = cqlx) onde & é a estimativa da palavra código 


transmitida e r é a palavra recebida. Como visto, isto pode ser 


realizado maximizando E c;p O Pp, onde 9,=P(r,|1)/P(rp]0) é a 
razão de verossimilhança. Novamente a dedução dos resultados é sim - 
plificada assumindo que a palavra toda nula foi transmitida. 

Conhecida a densidade espectral de potência unilateral ! 
do ruído no2+ as componentes tnp) do vetor de ruído n são variã 
veis aleatórias iid gaussianas, np O Nto,ng/2). Novamente será utili 
zada a relação sinal/ruído do canal y = E/ngr ea relação sinal/rui- 
do por bit de informação transmitido hi ba E,/ng = Y/R. 

Das considerações acima segue-se que a probabilidade de 


erro por bloco & 
Punddo * P[e | c=0] (3-53) 


Dado que a ênupla toda zero “c 
verã erro na decodificação se e somente se 


alguma palavra côdigo e, i = DR 


Como os sinais transmitidos são pulsos de amplitude JE 


volts ou O volts, e o canal é gaussiano, segue-se que 


tim 1vET? 
P(rp|i) = exp (= Smmemee— i=0,1 (3-54) 
Ing Ng 
e portanto, 
n=1 n=1 
E c;p nd, = E z gp - to (3-55) 
2=0 9/2. 2m0 * 
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Como é suposto que a palavra toda zero foi a palavra '! 


transmitida, o vetor recebido *r = c+n = (rgtlyrec err) tem |! 


sua t-êsima componente dada por rp = np. 


Se são definidos eventos Apis MPR 
n-1 n=1 

A, Ôtie nl 1 e) É (3-56) 
2=0 t=0 * 


então & possível reescrever (3-53) como 


= PLU A) (3-57) 


se 
indo a P(A) (3-58) 
i=1 
n-1 
Por outro lado, L c.p = wí(c.), onde Wwí(c.) represen- 
£=0 1 ee * fast N 
á à Med 
ta o peso da palavra c,. A variável aleatória X, = L ec,n e 
med, pi t=0 de 


uma combinação linear de variáveis aleatórias gaussianas independen 
tes, o que acarreta em X, tambêm ser uma variável gaussiana | 9|, 


n 
com distribuição de probabilidade Ro Nto, w(c)-S. Portanto 2 


+00 
PAD = P(X, > 15 w(c,) = : EMA) ax (3-59) 
É ua W(c;) 
ou seja, 
BA ds QUE W(c,)/2n9) 
= o(/RY, W(c,)/2) (3-60) 


= 


A cota (3-58) é apertada para alta relação sinal/ruído'! 


[19], de modo que é possível escrever 
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e) 
tm 
H 


( 


"loca!" e P(A;) (3-61) 


AB 
ci i=l 


Contudo, somente palavras código com peso d. contribu- 
em substancialmente no somatório em (3-61), visto que para valores! 


suficientemente grandes de x, Q(x)=z expt=x"/2]; Logo, o com -— 


x 


portamento assintótico para a probabilidade de se decodificar incor 
retamente uma palavra quando a regra ôtima é utilizada em um canal! 


gaussiano com alta relação sinal/ruído ê expresso por 


AB pro lPpioco! * N(d) Q(y/RdY,/2) 

Considerando pulsos polares de amplitude +yE volts e 
-/E volts, o resultado obtido coincide exatamente com o comportamen 
to assintôótico determinado para o algoritmo de Hartmann-Rudolph na 


seção 3.1.4. 


CAPÍTULO IV 


SIMULAÇÃO EM COMPUTADOR 


Neste capítulo são apresentados os resultados de simu 
lações realizadas em computador digital (DEC-10 SYSTEM) para análise 
do desempenho dos dois algoritmos descritos no capítulo anterior. A 
necessidade desta simulação advêm da complexidade* do procedimento ! 
analítico da variação da probabilidade de erro (por símbolo ou bloco) 
em função da relação sinal/ruido, e da degradação que resulta quando 
decisão suave é utilizada quantizando em Q regiões os símbolos re- 
cebidos do canal. O desempenho de vários côdigos ê analisado, e em 
particular, os códigos de Hamming (7,4,3), (15,11,3) e (31,26,3) são 
considerados. O intuito & fazer um levantamento das curvas de proba- 
bilidade de erro versus relação simal/ruído para cada côdigo de blo- 
co analisado, em um dado canal sem memória. Isto torna possível, por 
exemplo, que seja feita uma visualização da melhoria provida pelo '! 
uso de decisão suave com relação a decisão abrupta e com relação ao 
número de regiões de quantização utilizado. 

' 


Os programas foram escritos na linguagem FORTRAN 10 


e se encontram listados no apêndice A. 


*Veja apêndice B. 
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4.1 - ANÁLISE DO DESEMPENHO DO ALGORÍTMO DE HARTMANN-RUDOLPH 


4.1.1 - Considerações Gerais 


O desempenho do algorítmo de Hartmann-Rudolph será ana- 


lisado para vários côdigos de bloco lineares, quando os sinais É 


transmitidos são pulsos liga-desliga sujeitos a ação de um ruído ! 
aditivo gaussiano. Em particular, resultados para os códigos de 
Hamming são apresentados. 

Na simulação, primeiro é especificado o código linear ! 
binário C(n,k,d) fornecendo sua matriz de verificação de paridade 
[H] , ou o polinômio gerador g(x) no caso de códigos cíclicos. A 
seguir é fixado o número de regiões de quantização a ser utilizado! 
na simulação, usualmente uma potência de dois. A análise é realiza- 
da considerando 2,4,8 ou atê 16 regiões de quantização. 

A simulação ê feita transmitindo-se a palavra toda zero 
em um canal com ruído aditivo gaussiano, gerado de acordo com o mê- 
todo polar |22|. Não hã nenhuma perda de generalidade em se assumir 
que a palavra transmitida foi a palavra toda nula |19|. 

Para os códigos de Hamming analisados, a opção de utili 
zar a palavra transmitida como sendo a palavra toda um, forneceu |! 
praticamente os mesmos resultados obtidos quando a palavra transmi- 
tida foi a toda nula, mostrando que o ruído gerado apresenta sime - 
tria. Os resultados para o côdigo (7,4,3) tambêm foram verificados" 
considerando os bits do bloco de mensagem gerados de acordo com uma 
sequência pseudo-aleatória |11). 

As probabilidades de erro por símbolo e por bloco são 
então. estimadas para cada valor de relação sinal/ruído (em decibe- 
is) transmitindo um número suficientemente grande de palavras côdi- 
go atravês do canal e utilizando a regra ôtima de decodificação no 


receptor. A estimativa ê feita utilizando a frequência relativa de 


. 
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ocorrência de erros na saída do decodificador. 
A utilização da versão fraca da lei dos grandes números 
de acordo com o teorema de BERNOULLI | 9| assegura que os resulta - 
dos obtidos convergem para os valores reais das probabilidades ,pois 


se n,/n é a frequência relativa da ocorrência de erros, então 


n 
á ] 1 
P = 2 el£ -—— Para todo e > 0. - 

| | n | e ” âne ea 


Deve ser notado que o número de palavras a serem trans- 
mitidas pela fonte deve ser aumentado à medida que a relação sinal/ 
ruído cresce. Tal aumento de relação sinal/ruído acarreta uma dimi- 
nuição no número de erros provocados pelo canal e consequentemente 
a interpretação frequentista utilizada para o cálculo das probabili 
dades de erro poderã não mais ser válida. Pode ser facilmente de - 
monstrado |33| que para uma variável aleatória com distribuição bi- 
nominal Q3 (n,p), um estimador não enviezado para p & a frequên - 
1 


cia relativa p =k/n, o qual resulta em um espalhamento relativo 


em torno do valor esperado da estimativa dado por 


O (4-2) 
E(p) VR p 


Se p representa a probabilidade de exro por bit na 


saída do decodificador, e q a probabilidade do evento complemen - 


tar, então para p<s< 1 segue-se q=l, de modo que 


a a (4-3) 
p|3/E1P)| 


Por exemplo, se a probabilidade de erro &ê da ordem de 


Emi a - : = i 
10 e & desejado um espalhamento relativo de 10%, entao a simula - 


= 5 
ção deve ser realizada para n da ordem de 10 bits transmitidos . 


Desta forma, para SNR = 4dB é suficiente que o número de bits R 


. 


ETD= 


transmitidos seja da ordem de mil bits, enquanto que para -SNR = 12 
dB é necessário que este número seja cerca de 100 mil bits, para as 
segurar a validade dos resultados encontrados. 

Como a palavra decodificada pelo algoritmo não ê neces- 
sariamente uma palavra código, surgiu a idéia de escolher como esti 
mativa da palavra transmitida a palavra código mais próxima em ter- 
mos de distância de Hamming da palavra decodificada. Contudo, os re 
sultados obtidos como o emprego de tal procedimento degradaram sen-= 
sivelmente a probabilidade de erro por símbolo, enquanto que a me + 
lhoria na probabilidade de erro por bloco foi praticamente nula. 

O tempo de CPU requerido para realização das simulações 
é razoavelmente grande, como seria de se esperar, não chegando con- 
tudo, a tornar-se proibitivo, exceto para SNR's muito elevadas. Con 
“tudo, para altas relações sinal/ruído, o comportamento assintótico ! 
do algoritmo é conhecido. 

De posse dos resultados obtidos, foram traçadas as cur- 


vas de probabilidade de erro versus relação sinal/ruído mostradas a 


seguir, nas páginas 108 a 124. 


4.1.2 - Curvas de Desempenho 


As curvas de desempenho para o algoritmo de Hartmann - 
Rudolph obtidas atravês de simulação em computador digital estão 
apresentadas nesta seção. É admitido um canal perturbado por um ruí 
do branco gaussiano aditivo em todos os casos. O número de regiões" 


de quantização é abreviado por NQ e PBIT denota a probabilidade 


de erro por símbolo, enquanto que PBLOCO denota a probabilidade de 


erro por palavra. 


o” 7 | 


1:E=1 


iso 


Figura 4.1 - SIMULAÇÃO PARA O CÓDIGO(7,4,3) 
Curvas PBIT x SNR 


1 - Decisão abrupta NQ = 2 


2 - Decisão suave NQ = & 
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Figura 4.2 - SIMULAÇÃO PARA O CÓDIGO (7,4,3) 
Curvas PBLOCO x SNR 
1 - Decisão abrupta NQ = 2 


2 - Decisão suave NQ = & 
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Figura 4.3 - SIMULAÇÃO PARA O CÓDIGO (7,4,3) 


Curvas PBIT x SNR 
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8 E 
E Figura 4.5 - SIMULAÇÃO PARA O CÓDIGO (15,11,3) 
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Figura 4.6 - SIMULAÇÃO PARA O CÓDIGO (ss 3 
Curvas PBLOCO x SNR 
1- NQ=2 
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SIMULAÇÃO PARA O CÓDIGO (31,26,3) 


Curvas PBIT x SNR 
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SIMULAÇÃO PARA NQ = 
Curvas PBIT x SNR 
Côdigo (7,4,3). 
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SIMULAÇÃO PARA NQ = o 
Curvas PBLOCO x SNR 
Código (7,4,3) 
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SIMULAÇÃO PARA NQ = 2 
Curvas PBLOCO x SNR 
Código (7,4,3) 
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SIMULAÇÃO PARA NQ = 4 
Curvas PBIT x SNR 
código (7,4,3) 

código (15,11,3) 


Código (31,26,3) 
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SIMULAÇÃO PARA NQ =. 


Curvas PBLOCO x SNR : 


Código (7,4,3) 


Código (15,11,3) 
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Figura 4.14 - SIMULAÇÃO PARA NQ = 8 =121= 
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SIMULAÇÃO PARA NQ = 
Curvas PBLOCO x SNR 
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4.2 - DESEMPENHO DA DECODIFICAÇÃO UTILIZANDO MÁXIMA PROBABILIDADE A 
POSTERIORI 


4.2.1 - Considerações Gerais 


Uma avaliação serã apresentada do desempenho do segundo 
algorítmo ótimo descrito no capítulo anterior (sec.3.2)para vários ! 
côdigos de bloco lineares, supondo sinais transmitidos em pulsos de 
RF liga-desliga atravês de um canal ruidoso. 

O algoritmo fará uso da treliça associada ao código pa- 
ra simplificar as operações de decodificação. Inicialmente um códi- 
go de bloco linear em GF(q) é especificado atravês de (n,k,d) e da 
sua matriz de verificação de paridade. A treliça associada ao códi- 
go é então gerada de acordo com um procedimento melhor esclarecido" 
na seção 4.2.2. 

As palavras código transmitidas pela fonte são geradas 
aleatoriamente e de tal forma que possuam uma distribuição de proba 
bilidade uniforme. Deste modo, a hipótese requerida para otimalida- 
de do algoritmo & obedecida. O ruído na representação de banda es - 


k Serão 


treita ê novamente simulado utilizando o método polar |22 
abordados apenas códigos lineares binários, muito embora a treliça" 
possa também ser gerada para códigos de bloco lineares multiníveis, 
As probabilidades de erro por palavra e por bit quando este algorit 


mo de decodificação é utilizado na recepção são determinadas para ! 


cada valor de relação sinal/ruído (em dB), transmitindo um número ' 
suficientemente grande de palavras código atravês do canal ruidoso. 
A estimativa utilizada para determinação destas probabilidades e 
novamente a frequência relativa da ocorrência de erros, de modo que 
as considerações sobre convergência dos resultados discutidas na se 
ção anterior ainda permanecem válidas. Com relação ao tempo de CPU 


requerido, considerações similares as enunciadas na seção anterior" 


tambêm são aplicáveis. 
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Qs24 do Implementação Computacional da Treliça 
Serã agora definida uma matriz associada à treliça gera 
da por um código de bloco linear em GF(qg), com o intuito de permi - 


tir uma implementação da treliça sob o ponto de vista computacional 


E 


Seja [M] uma matriz com q trt) linhas e n colunas, 


onde cada elemento desta matriz ê um vetor de q posições, ie, 


M,1 m,,2 cemererco My 
Ro al LR) neces. om 
[M] = Ê (4-4) 
m En edo O RÁ : m n-k 
= n-k = n-k a 
el q 12 q sn 
EXis qu 
com m, RI aa ae ) 
=i,k QE "“g-1 Ll<k<n 
-k 


Ep E io E 


Portanto, o número de posições de memória necessário pa 
ra definir uma treliça de um código linear (n,k,d) em GF(g) será 
de 


n Ea =n qe posições. (4-5) 


Cada vetor coluna de [M] estã associado a uma profundi- 
dade k > 0 da treliça, e cada elemento Mi deste vetor coluna ! 
estã associado a um nô da mesma. Cada elemento my Por sua vez , 


possui sua coordenada u, igual a numeração correspondente ao nô 


k da profundidade anterior tal que sp (k) = s; 0-1) (4 hp. Deste 


* 
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modo, cada vetor ua da profundidade k indica quais os nôs e ele 


mentos Gs e GF(q) que deram origem ao nô ao qual o vetor Ma es 
tã associado. Fixados um dado nô NO e uma profundidade PROF, ! 
MATRIX(NO, PROF, ALFA), ALFA = 0,1,....,q-l define uma q-upla cor 
' 


respondente a m A primeira posição da gq-upla significa 


NO, PROF * 


que deve ser considerado que o nô em questão foi atingido atravês ' 
de a = 0, enquanto a segunda posição significa que o ramo considera 
do se dirige para o nô em questão atravês de a = 1, e assim por di- 
ante. Os valores indicados na q-upla fornecem os nôs na profundi- 
dade anterior que se ligam com o nó em questão. É assumido que (o) 
valor zero indica que nenhuma ligação deve ser considerada. Como ' 
exemplo, o valor de MATRIX(3,2,0) fornece qual o nô na profundidade 
1, o qual atravês de q = 0 atinge o nô 3 na profundidade 2, 

O procedimento utilizado para expurgar a treliça asso - 
ciada ao côdigo de bloco considera inicialmente a construção de du- 
as treliças auxiliares, A primeira delas, TRELLI 1, é gerada a par- 
tir da profundidade k = O atê a profundidade k = n; enquanto que 
a segunda, TRELLI 2, ê gerada iniciando na profundidade k =n atê 
a profundidade k = 0. Os caminhos que coincidirem nas duas treli - 
ças fazem parte da treliça expurgada, em caso contrário o caminho ê 
apagado. Este procedimento & melhor compreendido quando apresentado 
nos exemplos mostrados a seguir. 

Exemplo 4.1 - Caso binário. 


Dado o código de Bloco (5,3,2] com matriz de verifica - 


ção de paridade [H] 


DE SDS RP 
[8] = 


Lt 7 o 


A matriz [M] que representa a treliça associada ao cô- 
digo é construída de acordo com os passos mostrados abaixo, 


Tem-se n=5, k=3, n-k=2, d=2, q=2 e qP-k do 


» 
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à 
Estados 


K0  K=l K=2  K=3 K=4  Ks5 TRELLI 1 


6,0) 6,9) 6,4 0,3 0,2 
0,1) 2,0) 2,3 2,4 2,1) 
0,9) [0,1] 8,7 B,) B,9 
0,0 [0,7 &,) [4,2 [4,3 


1= (8) 
2= (9) 
3= (1) 
4= (1) 


TRELLI 2 


6,3 0,3 0,9 0,3 0,2 
2,7 2,9 2,3 2,9 [,9 
6.9 6,] 6,3 D,9 [0,9 
6,3 6,9 [4,7 ,9 D,9 


K=0 K=1 K=2 Kk=3 K=4 K=5 


Finalmente a treliça expurgada associada ao código binã 


rio (5,3,2) é encontrada usando TRELLI 1 e TRELLI 2. 


TRELLI 


0,9) 6,9 0,9 0,3 0,2 
0,7) 2,9 [2,3 2,9 [0,9 
0,9 [0,7 [2,2 [0,9 [0,9 
0,9 [0,27 [4,1] [0,0 [0,9 


Exemplo 4.2 - Caso multinível. 


Como um segundo exemplo serã considerado um côdigo de 
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bloco (4,2,3) definido em GF(3), com matriz de verificação de pari- 


dade dada por 


Os parâmetros são: 


k 


n=4, k=2, n-k=2, d=3, q=3 e es =9, 


A treliça associada a este código linear e a matriz [M] 


representando a treliça são mostradas: 


1= (8) 
2= (9) 
3= (9) 
4= (0) 
5= (1) 
6= (3 
7= (5) 
8= (1) 


= (3) 


K=0 


Profundidade 
Estados 


k=1 


(1,0,0) (1,0,9 [1,7,4] (1,3,2 
[0,0,0) [0,0,8] [2,8,5) [0,0,0] 
[0,0,0] [0,6,0] [3,9,8) [0,0,0) 
[0,0,0] [0,0,1] [0,0,0] [0,0,0) 
[0,0,9 [0,8,9) [0,0,0) [0,0,9 
[0,1,0)] [6,0,0] [0,0,0] [0,0,0] 
[0,0,0] [0,1,0) [0,0,9 [0,0,9 
[0,0,] [8,0,0) [0,0,0] [0,0,9] 
[0,0,9 [0,0,8] [0,0,0] [5,0,9] 
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A geração da treliça permite a implementação de diver - 
sos algorítmos para decodificação de códigos de bloco lineares, al- 


guns apresentados nesta tese. 


4.2.3 - Curvas de Desempenho 


As curvas de desempenho para o algoritmo de maximização 
da probabilidade a posteriori das palavras código obtidas através ' 


de simulação em computador digital estão apresentadas nesta seção . 


As considerações descritas na página 107 sobre a notação utilizada" 
novamente se aplicam. Estas curvas permitem uma visualização da me- 
lhoria obtida com uso de decisão suave. O caso onde decisão suave é 
utilizada quantizando em Q regiões as amostras recebidas do canal 


é também considerado. 


E) 


Figura 4.18 - SIMULAÇÃO PARA O CÓDIGO (7,4,3) =131- 


RU =4) 


Curvas PBIT x SNR 


i 


1 
[8 


1 - Decisão abrupta NQ 


| 
8 


2 - Decisão suave  NQ 


8 


& 


“Ba 


1 


cgi É 


1 


dos yo SS VLS 


& 
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Figura 4.20 - SIMULAÇÃO PARA O CÓDIGO (7,4,3) Ao 
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Figura 4.21 - SIMULAÇÃO PARA O CÓDIGO (7,4,3) 


Curvas PBLOCO x SNR (Deteção por envoltóôria) 
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Figura 4.23 - SIMULAÇÃO PARA O CÓDIGO (7,4,3), 
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CAPÍTULO W 


CONCLUSÕES 


5.1 - ANÁLISE DOS RESULTADOS 


O intuito principal deste trabalho & o de apresentar 
e analisar procedimentos de decodificação de códigos corretores de 
erros, em função da melhoria proporcionada no desempenho de sistemas 
codificados. Como visto, tal melhoria é obtida através do uso da in- 
formação probabilística associada as amostras recebidas, e em geral" 
& conseguida as custas de um aumento na complexidade do sistema. 

No capítulo II, após a revisão sobre códigos linea - 
res, foram introduzidos procedimentos subôtimos que empregam têcni - 
cas de decisão suave (sec.2.2), os quais proporcionam uma melhoria ' 
considerável no desempenho do sistema, como mostrado nas curvas cor- 
respondentes as figuras 2.4, 2.6 e 2.8. Eles são muito atrativos le- 
vando-se em consideração que o uso de regiões de quantização & ade - 
quado por facilitar a implementação utilizando técnicas digitais. O 
procedimento proposto por Wolfenson-Rocha & de grande importância |, 
visto que é um dos poucos esquemas de decodificação que assume ser 
, 


desconhecida a distribuição de probabilidade da fonte. A melhoria 


provida pelo uso desta têcnica encontra-se analisada com detalhes na 


referência |7|. O desempenho do algoritmo de decodificação por dis - 


. 


SB 


tância generalizada II foi analisado atravês de simulação em compu- 
tador, e comparado com o desempenho do sistema empregando decisão ! 
abrupta, resultando nas curvas das figuras 5.1 e 5.2 mostradas a se 
quir. 

Para os procedimentos ôtimos descritos no capítulo 
III, o estabelecimento da regra de decodificação permitiu a determi 
nação do comportamento assintótico para alta relação sinal/ruído.As 
simulações realizadas para estes algoritmos estão apresentadas nas 
seções 4.1.2 e 4.2.3 e a obtenção destas curvas são de importância! 
fundamental neste trabalho. A partir de cerca de 12dB, a aproxima - 
ção utilizando a expressão assintôtica fornece excelentes resulta - 
dos, e & acima deste valor que a simulação se torna praticamente in 
viável. Deve ser observado que jã existiram problemas na convergên- 
cia da estimação da probabilidade de erro quando utilizando o códi- 
go (31,26,3), isto porque a regra de decodificação descrita pela 
equação (3-2) torna-se excessivamente complexa* aumentando de modo 
sensível o tempo necessário para decodificar um bit. Este fato estã 
refletido em algumas curvas da seção 4.1.2 para os maiores valores" 
da relação sinal/ruído. 

Quando o número de níveis de quantização ê aumenta 
do, obtêm-se informação probabilistica mais detalhada para ser uti- 
lizada pelo decodificador. Entretanto, um aumento no número de ni - 
veis de quantização significa um crescimento na complexidade dos 
circuitos detetores. Em todos os procedimentos descritos, observa - 
se que a medida que o valor de NQ aumenta, o ganho incremental na 
relação sinal/ruído para se obter uma dada probabilidade de erro & 
cada vez menor, de modo que grande parte da degradação que ocorre ! 
quando NQ = 2 pode ser superada sem que a complexidade do sistema 
se torne proibitiva. A informação adicional obtida por quantizar as 
amostras recebidas em mais que 16 níveis acrescenta muito pouco d 


âquela obtida para 8 níveis. 


* Compare este caso com o descrito no exemplo da página 83. 
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Figura 5.1 - SIMULAÇÃO PARA O CÓDIGO (7,4,3) 


Curvas PBIT x SNR 


1 - Decisão abrupta NQ = 2 


2 - Decodificação por distância generalizada II 
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“ Figura 5.2 - SIMULAÇÃO PARA O CÓDIGO (7,4,3) -140- 
Curvas PBLOCO x SNR 


1 - Decisão abrupta NQ = 2 


2 - Decodificação por distância generalizada II 
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As curvas citadas na página 105, usadas para veri- 
ficar a simetria do ruído Gaussiano gerado são apresentadas a se - 


guir. 
5.2 - COMENTÁRIOS 


Estas técnicas descritas são de particular interes 
se para sistemas de comunicação via satélite, em comunicação mili - 
tar e em comunicação espacial de um modo geral. Uma importante e 
atrativa aplicação estã no emprego da decodificação probabilística" 
em sistemas de comunicação utilizados para transmissão de dados em 
canais de HF, de modo a assegurar uma boa confiabilidade à transmis 
são. 

Neste aspecto, ê sugerida uma continuação deste '! 
trabalho, propondo-se análise e implementação dos algoritmos estuda 
dos, aplicados para canais de HF. 

Com relação ao procedimento proposto por J.Wolf ' 
descrito no capítulo II, deve ser mencionado que & um procedimento" 
ótimo, equivalente ao algoritmo de maximização da probabilidade a 
posteriori descrito no capítulo III. Contudo, a formulação matemáti 
ca no último caso, alêm de maior rigor, permite a interpretação de 
"correlação generalizada" e o estabelecimento de cota para a proba- 
bilidade de erro por bloco. O mais importante ê que permite o uso 
de decisão suave usando NQ = 22 regiões de quantização, o que de 
acordo com a simulação resulta em considerável melhoria com relação 
a decisão abrupta e uma simplificação enorme no decodificador. Ade- 
mais, comparando-se as expressões utilizadas na decodificação, vê - 
se que a segunda apresenta uma vantagem computacional. 

No caso de côdigos de bloco, fica claro que o algo 
rítmo de Hartmann-Rudolph pode ser implementado utilizando a treli- 


ça associada ao código dual. A implementação computacional da treli 
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ça na forma descrita na seção 4.2.2 e nas subrotinas das listagens do 
apêndice A, permite a simulação de diversos tipos de decodificadores 
de códigos de bloco. 

Algumas versões dos programas estão apresentadas 
nas listagens do apêndice A. Entretanto, devido ao tempo de CPU re:- 
querido, os programas não podem ser codificados para um caso geral , 
de modo que devem ser processadas algumas pequenas alterações para 
cada caso específico. 

Ainda como sugestões para pesquisas futuras, de- 
ver-se-ia considerar o uso de microprocessadores para implementação! 
dos decodificadores Ótimos propostos no capítulo III, permitindo in- 
clusive o emprego de subtreliças. Com relação ao algoritmo de decodi 
ficação por distância generalizada II, seria interessante tentar de- 
senvolver uma cota para a probabilidade de erro. 

É interessante observar que a decodificação por 
distância generalizada II resulta em uma palavra onde a mêdida dos 
coeficientes de confiabilidades condicionais dos dígitos & máxima . 
Este fato encontra-se demonstrado no apêndice C. 

A medida de confiabilidade associada a uma clas- 
se pode ser determinada mesmo no caso onde a partição do espaço de 
observações resulta em classes de confiabilidade não simétricas, e 
pode ser calculada no caso multinível. É possível conjecturar que 
exista uma maneira Ótima de combinar os coeficientes de confiabilida 
de p 6) de modo a obter um algoritmo que minimize a probabilidade" 
de erro por símbolo (similar ao descrito no capítulo III) que possa" 
ser aplicável a códigos de baixa eficiência e que seja simples no 
caso multinível. Uma posterior extensão para decodificação de côdiv- 
gos convolucionais seria bastante útil, pois neste caso os códigos ' 
utilizados geralmente são de taxa 1/n. 

As têcnicas de decodificação apresentadas nesta" 


tese se mostram muito poderosas e práticas em canais que tenham com- 


. 


mL4S- 


portamento aproximadamente semelhantes a canais sem memória. Com o 
custo do hardware digital diminuindo rapidamente, o uso destas têc- 
nicas crescerá largamente em aplicações práticas em situações de 


controle de erros. 


APENDICE A - 146 - 
c UF PEL/ULS= CUMENE 
(o DECODIFICACAU DE CULIGUS LINLAKES 
c ALGURITSU DE hARTMANAN É HyUDbULEN 


INteEGEI ALL DIL,BOLUND,ChECR,D, LMLN,G, hn, DEZ, TA 

DIMENSLON DOUND(4, 4) V(3IZ,31),G(IZ),HC3Z),PICIO), PUCIO),K(C16),KUCS3 
=1),X132) 

CUsmuN DoNCobhs Lou 

DATA 6/3240/,X/32*%0/ ,VES/"YES!Z 

EXF(X)=1 -ERECILX) 

SuZz5UkRIT(2.) 


(6 ENTRALA DE DALOS à FHLUCA CUUt 
nnxitetD,11) 

11 FUKRATCIHO, "ENTER norqL=",/) 
B=20 


UPES LUnLI=20,FiLE='HARD DAT!) 
READ(M, +)NC, RC, Dula 
NENC-RC 
IF(NC,GT.32.UK NGT.S)STUP 
n=2++ 
ACPI=AC+1 
UCPI="NC+1 
nPl=0+1 
nLizivel 
nLI=“=-3 
nKIIE(S, 22) 
.! FORMATLIHU, !CICLE CUDÊE?Z YES Uk nO!) 
nRrAvLis,3I)TES de 
33 FORMATÇAS3) r 
AFMTESTE heçitS)EU TO 1000 
kUTISA PZ CUDIGUS CÍCLICOS, LPUOLINOMIO GeKAVUR) 
(1,4) => (3,140) 1,0,1,1 
(15,11) => (4,140) 1,0,0,1,1 
(31,26) E E542,0) 1,0, V,1l,UVçl 
vRITE(S,44) 
44 FURSATCIHO, "ENTER G(A)=",/) 
KEAD(M, F)CGCI), 1=),8PI1) 
A(t1)=1 
A(uCP1)=1 
LDU Z00 n=1,nCra 
IF(X(A) EV. U)60O TU 100 
HtK)=1 
CALL SMODZCA,G) 
CALL SHJFICG) 
GU FU 200 
10U CALL SHIFTCG) 
E 2UU CuUniinUE 
bu 3uuv 1=1,34 
3VU CHECh=CHECK+ACÇI) 
1+(ChHeCReNE U)STUP 
LU 500 j=1,“Ll 
DU 4UU J=1, 40 
4UU DCI, NCPI-J)=Nty) 
CALL sSHIF ICH) 
SUU CusTisuUr 
GU TU ZUUu 


oqnoQn 


E kOLINA PZ CJLIGUS DE BLUCU NAN=-CICLICOUS, 
JUUU anITE(S,55) 
o FURBALÇIHO, fEntek U(n4)!,/) 

/775 DU oUU 1=1,5 
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bo FURMAT(3IZCIL, AD) 
oLU CunIVInUL 
E GERANDO PALAVKAS DU CULDIGU DUAL. 


20U0 h=u 
DO 700 1=1,wLl 
BUUNDEIL, I)=K+1 
141=]+1 
DU 700 J=IPI,N 
K=K+1 
CALL SCham 

190 CunTisvL 
Jr (ineLI.3)6U LO 500 
lJenb=aL1 
DU cOU L=2,NLI 
lEnDb=1END-1 
DO vOUO 1=1,1t5D 
BUUND(L, 1)=n+1 
JU=AUUNDIL=1,1+13 
JE=BUUNDIL=-1 NLI) 
DO GOV J=JU Ut 
K=n41 
CALL, SCRAM 

B0U Cuntinue 


E PNSIASSEFESRALDAFARASIRASLESESLSENSFRRAALASAFALERALEAPRRARR AREA 
I=(021n-1)/2 
A=1U 
ALL=O 

cs nUmtku LE kKEGÍIUVES Dt GUANTIZACAU = nQ 


DO 9U JhN=2,4,1 
NO=2*+iiy 
mRITE(3,1])ivd 
1 FURBAT(IHI,49X, URIVERSIDADE FEDERAL DE PLkNAaMBUCU!,/,51X, 'MESTRAD 
-U Ex cnGtnHanIa ELETRICA!,/,30X,'4% DECODIFACACAU DE CODIGUS LINEA 
“RES - ALGURILIMU DE HARTMANN & KUDOLPH F51,/,5X, !hu=!, 14) 
nol=no-1 
CA VAKIANDO SIGNAL-TU-NUISE FATIU, 
DU 50 SNk=11c,iZeção 
NX1=0 
2=10,**(Stik/20.) 
E K=1 (LOW NOISE) 
B=Z*4(»U41)/(NG+SUZ) 
PI(1)=0.S*ERFCCS) 
B=4/(unGtSU4Z) 
PULI)=0.5*L1,.4ERECb)) 
E K> nG-1 0 (HIGH nUISt) 
PI(ny)=PULI) 
PULNG)=PI(I) 
c I<k<no=1 
1 (nuçto,Z2)6U TU 30 
DU 30 n=2,NuÚl 
BI=Z+(NQ=N+1)/(uutouz) 
B2=LE(Ny=-k)/(nOt5uZ) 
PiCN)=0.5*(LHE(BI)-ERECEZ)) 
PUCNG=-n+41)=PÍCA) 
30 CONTINUE 
é) CALCULU DOS RU!'S 
DU bu 1=i Nú 
RCI)D=C1,=PICI)/POLI))/C1 4PICID/FUCI)) 
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THE=U. 
1Pr=u. 
E» 
B1LT=(10*%4 )$NC 
STU=0. 
x SiIMULAÇAU DUO ALGURITMU 
DU 40 KDIG=1,EIT,NC 
NYZ=u 
DU 10 J=1,nC 
G GERANDO KULDU GAUSSIANUS 
15 vlz2*KAN(5U)-1 
V2=2*kKAN(34)-1 
S=VI+AVI+VZHVZ 
LIFE (5.GE,1) GO TU 15 
V=VIASURIC-ZHALUG(S)/5) 
I=SIGmA*V 
c POSSIBILIDADES: 
IF (reGL.A/2. )nY2Z2=Nny241 
c IFCY.LEÇ-A/2. )NYZ=Ny2+1 
E INTERVALO UNDE CAIU À AMOSTRA =hY 
NI=LALLFA+Y)ENU/A+1 
E rROB Y DIVISiBLE BY A/NQ. IS ZERO 
IFCnT.LT.I)NY=1 
LE CNY GT NQI)NT=ZNQ 
RUCI)=RENY) 
10 CUNTINUL 
nxl=nNxI+NY2 
IF (NyZeLE SI) GO TG. 40 
DU 35 1=1,nC 
SP=0. 
DO 25 J=1,w 
PROD=1 
DU 20 L=1,nC 
DELTIL=1 
R JF CG NE.L)DELTILSO, 
IFCD(J,L).NESDELTIL)PROD=PRODFRUCL) 
20 CONTINUE 
25 SP=PROD+SP 


C DECISAU : SP>O, ESTIMAR Cº=0, SP<O, ESTIMAR 


(6 POSSIBILIDADES: 
IF(SP.LE.O)TBE=TBE+1 
Cc JFCSP.GT.O)TBE=TBE+I1 
35 CONTINUE 
IF(TBE.GI.STO)TPE=TPE+] 
STD=VEF 
40 CONTINUE 
c CALCULANDO PROBABIÍ. IDADES DE ERRO, 
PEB=Tot/BIT 
PEP=TPEANC/BIT 
E SAIDA DOS DADOS 
WRITE(3,2)SNkR, THE, PER, TPE,PEP, NY 


= AB 


2 FORMAT(IHO,15X, !SNR=",E14,8,4X, "TBL=",E14,8,4X, "PEBS!,E14,8,4X,!T 
-PE=1,E14.8,4X, "PEP=",614,8,/,10",37X, !NYI=",16) 


50 CONTINUE 
STOP 
END 
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SUBROTIKA P/ CALCULO DE EkF - APROX, POLINOXITAL 
FUNCAO DISTRTBNTCAD SORKAL DE PROBABILIDADE dado 
ERRO /E(X)/ > 7.5 E-08 
vEES - DEPARTAPENTO ENGENHARIA ELETRONICA E SISTEHAS 
PuULCciION ERFCC(Z) 
X=ABS(SGET(Z.)%Z) 
SGET(1/(2%3.141593)) 

2316419 

b1= 3193615 

- b2=-, 3565638 
83=1,781476 
B4=-1.821256 
E5=1,3302740 


Qonnn 


Dc 


T=1/(1+P*X) 

ERFC= QFEXP((=-x**2)/2) 
ERFC= ERFCATx(RISTE(B2+TH(B3+TH(B4+T*HEB5)))) 
ERFC=2,*HERFC 
RETURE 
END 
SUBRDUTINE Sz0p2 (X,6G) 
INTEGER X(32),G(32) 
DO 1 1=1,32 

1 XCI)=IABSC(XCT)-GCI)) 
RETUEN 
EuD 
SUBEROUTIKRE SHIFTC(S) 
IuTEGER G(32),ST0l,STUZ 
STOI=G(1) 
b6 1 1=1,31 
STU2Z=G(1+1) 
G(1+1)=ST01 
STOi=STO2 

1 CONTINUE 
G(1)=ST01 
RETURS 
END 

(a SUBROUTINE 

SUBROUTINE SCRAN 
InTEGER D 
DIZENSIUN D(32,31) 
COFBON D,NC,K,Touy 
DO 1 L=1,NC 

1 D(K,L)=ABS(D(I,L)I=D(J,L)) 
RETUPE 
ED 
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UFPE/DES - CoMENF 
DECUODIFICACAN DE CODIGOS LINEARES 
ALGuRITAO DE DECODIFICACAO POR MAXIMA PROBABILIDADE A POSLERIORI 
COrRMON GE,WN.NPS,J 
INZEGER C(7),GF.GFO,HAT(7).HC(3,7) 
DIKENSIOw MATRIX(8,7,02:1),MATRIZ(8,7,0:1),PETRIXC(B,7,0:1) 
DIKRENSION R(7).STAR(7),PI(16),PO(16),ALCGFICI6) 
ERE(A)=1.-ERFCCX) 
S02=S0RT(2.) 
nQzie;NQI=NQ-1 
N jS=) 0r%4 
“RITE(5,33)N,K.GFO 
“RITE(5,11) 
OPEN(UNIT=20,FILE=!H.DAT!) 
DO 1 JJ=1,9 
READ(20. 7) CHCJJ, RN), 
«RITE(S, 10) CHOUJ NS) NNZ 
CALIL TRELLICh,MATRIX, MATRIZ) 
CALL EXPUR(MATRIX. MATRIZ) 
“RITE(5S,22) 
DO 2 K1=1,NPS 
“RITE(S,100)(C(UATRIXC(K1,K2,K3),K3=0,6F),K2=1,N) 
DO 5 SNRS0.0i2çwls 
TDr=0. 
TPE=0. 
2=10.**(SNR/20,) 
B=7*N01/(NQ+SG92) 
P1(1)=0,5*ERFC(B) 
2=2/(NQ*509) A 
PO(1)=0.5*(1.+FRF(CE)) 
PICNQ)=POC1) 
PO(NQ)=P1C1) 
1F(NQ.LE.2ISTUP 
DO 17 IREG=2,N01 
B1=Z*+(NO-IREG+1)/(N9%S0Z) 
BZ2=Z*(NO-IREG)/(N0*502) 
PI(IREG)=0,5BX(ERF(BI)-ERF(B2)) 
PU(NO-IREG+1)=PICIKEG) 
7 CouNTINUE 
DO 231 IREG=1.N0 
31 ALOGFI(CIREG)=ALOG(PICIREG)/POCIREG)) 
LO 4 NUM=1,N4nRDS 
CALL “ORDS(MATRTX,C) 
XSNR=SuR 
CALL NUTSE(XSNRE,NoR,Co STAR) 
DO 24 KX=1,N 
E) STAR(KXA)=ALOGFI(IFIXCSTARCKX))) 
CALL CURREL(“ALRIX, PETRIX, STAR) 
CALL SURVIV(rATRIX,PoTRIX.HAT) 
JIERRURDO 
DO 3 KX=14N 
IFREGR=LERROKR+IR=SCC(KYXI)-HAT(KX)) 


lp TEL=IBE+IER 
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44 


JECTERROR, NE.O)TPE=STPRES41 


COuTINUE 


PED= 


TRL/(NABANkDS) 


PEFSTPE/iIVORDS 
“RITE(S,44)SwR, PER, PEP 


U 


Form 
F Orr” 
FOR! 
FOnt 
FOR: 
-DO EM 


-ARES 


MiTINUE 


AT(1H .8X.10(11,1X)) 

ATCLHO SD dO CI LITE A 1 TT 1 2X 

AT(1hnO,BX, !'KATRIZ DE VERIFICACAO DE PARIDADE H!) 
AT(1h0,//,6X, !MATRIZ ASSOCIADA A TRELICA F4*%1) 

AT(1h1,44X; 'UNIVERSIDADE FEDERAL DE PERNAMBUCO!,/,47X, ! MESTRA 
ENGENHARIA ELETEICA!,/,30X,!** DECCDIFICACAO DE CODIGOS LINE 
- ALGORTIMO DE POCHA JR 451,///,9X,!'CNDIGO DE BLOCO (!,12,'!,! 


“4120 1)"3K MSDBRE GEC!) 1)çÃ//2) 


FOR! 
STOP 


ERD 


AT(I1HO,5X,!SwR=1,F4.1,5X,'PEB=",E14,6,5X, !PEF=1,E14,8) 
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taste 


SUSKOTINA P/ CALCULO Dk EKE - APKDOA. POLIND4IAL Se 
FUNCAU DISILKIBUTCAO NURNAL DE PROBABILIDADE 

EREUD ZECA)/ > 7,5 t-08 

DELES - DEPARTAMENTO ENGENHARIA ELETRONICA E SISTEMAS 
FunCTIOA ERFCCZ) 

X=hES(SORT(Z,)%2) 

U= S967(1/(2+3,141593)) 

p= .2316419 

Bi= .3193515 

E2=-, 3565636 

B3=1.751478 

B4=-1,821256 

D5ô5=1.330274u 


“aooaa 


T=1/(1+P+X) 
ERFC= OFEXP((-1A**42)/2) 
ERFC= ERECATICBI4TA(B24TIÇEI+TAÇES4)2B5)))) 
ERFC=2, fERFC 
RETURN 
ESD 
CHITRELLIEZNOCALESTRANDDRITKANDDEDEXPUKGONURDSA INODISEFECORRELSISURVIV 
E 
CH GtERanDO A IkKELICA NAG=EXPURGADA 
c 
SUBROGUTILNE IReLLIChA, MATEIX MATKIZ) 
CuUrrON GE, N,tPS,U 
InTeGER ALFA GE, CJ, A) EATERIACRES, à ,0:GE) MATRIZCNPS, às, O!GF) 
[cs GERasbLU à TRELLIS 1 (Em mATRIX) 
Cali, NUCALCI, |, d SBTRIA ÇÃO) 
LO 1 Kk=2,8,1 
BRKA=KAX 
Ni=kx-1 
bo 1 NAT=1 "PS 
mSAT="AT 
Istn=0 
DO 13 ALha=U, Gt 
13 ISUZM=ISUM+AMATRIXCANAL,FY,ALÉA) 
IFCISUM.EQ.0)GU TO 1 
CaLL NUCALCI ERA ENAT,MATRIX,H) 
1 CONTINUE 
C GERANDO À TRELLIS 2 CEM maTRIzZ) 
CALL NOUCALC2, 8, 1, SATRIZ,H) 
DO 2 kX=n-=1,1,-1 
MKAZRA 
RY=NX+41 
DO 2 mAT=1, Pê 
mNAT=5AT 
E DU Z ALFAZU,Gt 
IELSAIKRIZCANAT, II, ALÉA). LOS U)GU TU 2 
CALL NOCALCZ ha, PATRIZCSNAT,NI,ALEA) MATRIZ ÇÃO) 


2 CunilwUt 
kE Tui 
[MD 
c . 
Cez CALCULANDO US NUS - PkEtNCHIFENTO DA MATRIZ DA TRELICA 
c 


SUBROUTINE NOCALCET,KX AT SATRÇÃO) 

COMADA Gk,n,nPS,J 

InFEGEkR AT,ALFA ETCS),CTOS),GE,GFO,MATR(VPS, nn, DIGE)  HCJ,N) 
1S1G=(=1)+**("T-1) 
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GFO=GE4] a 


CALL TKRANDDBCGEFO, J, ATi) 

VU 2 abra=0,Gt 

DO 1 Ji=1,U 

LL=BT(UI)+ISIGHALFAFHCJI,RA) 

CrtuIi)="UDCLL,GFO) 

JFCCTCJI)D.LI O)CTCIID=CTCI1)A+AGEO 
1 CUusTINUE p 

CALL TRANBDCGEFO, Je NO,CT) 

JEÉCNT.EG.I)RATR(NO,KA,ALEA)=A1 

IF (NTLEG.Z)MATEÇAT,KA, ALEA)=NU 


2 CosTINUE 
ReTuRo 
ED 
É 
Cas IRANSFURSACÇÃO dl 
c 


SUBROUTINE TkAnBDCGFO,J, Db, BE) 
InTEGER LD, BB(D), GEO 

DD=0 

DO 1 Z=1,9 
DD=VD+BECJ+1-M)HGFO+*(M=-1) 


E CONTINUE 
Lo=vD+1 
RETURN 
Chu 
Ss 
c:4 TKANSFOKMACAUL Des 
c 


SUEROUTINt IKANDE(GEFG,J, BD, DE) 
INTEGER DD, BE(D),GEQ,4U 

u=PD=1 

DO 1 8=1,0) 
Bu(J+i=-2)=2UD(u,GFO) 


9=0/GFQ 
1 CONTINUE 
RETURN 
EnD 
E 
Ccr5 GERANDO A IRELICA EXPURGADA 
E: 


SUBROUTINE EAPURCSATRIX MATRIZ) 
COHMON GEN NPS,J 
INTEGER ALFA, GE LATKIX(NPS,N  USGE) MATEIZUNPS, N  DIGE) 
DO 1 NAJ=I,NPS 
DO 1 Kkx=1,N 
DO à ALFA=0,GF 
= JECMATRIACRAL, KA ELFA) it BATKIZCNAT, RA ALEA)) 
MATKIX (NAT, KX, ALF A)=O 
4ATRIZ(NAL,KX  ALFA)SO 


1 COxFINUE 
RETURN 
E END 
c ç 
Cre GSEkKANDU PALAVKAS=CULIGO ALERIULRTANENTE 
c 


SUKRKOLTIinE ALRDS(EATHIX,C) 
CussOiN GÉ,N,DPS,J 
IxIEGER ALÉA,CCN),GE,MALKIAXCUPS,N,02GÊ) 
Gru=ELOAMTUGSE4I) 
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Cru 


E) 
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nú=1 

DO 1 KX=n,1,-1 

ALFAZSFANCSO)AGEU 

JELALER EucGEU)STOP 
JECRALKIACNU, RX ALFA) EO. V)GU TU dO 
CCKX)=ALFA 

NOZMATRIXCNO, KA ALEA) 

Con) InuÊ 

RETURN 

Esb 


SIMULANDO CArRAL FOIDOSO - ENCONTRANDO VETOR RECEBIDO 


SUBROUTINE NULSECSWR,n,K,C,STAL) 
InTeGER Clin) 

DIMENSION K(b), STAR) 
nNo=16;NUl=nG=1 

h=1U. 

ZL=10.**(5S5k/20,) 
S)Graza/% 
VAR=SIGRAREZ. 

DO 2 1=1,u 
K=2.*RAN(S50)-1, 

Y=2,. *RAN(349)-1, 
S=ATX+Y*Y 
IF(S.GE.1.)GO tU 4 
V=5IGMALSORIC=-Z2,*ALUGCS)/5) 
AZA+V 

Y=1+v 

KCI)=CCI)*HA+A 
NISKCI)ANa/ A+] 

1F CRY LT.I)NTI 
JECNYSGTENQI)NI=NQ 
STAKCI)=NI 

CONTINUE 

RETURN 

Ent 


DETEEMINANDU À CURKRELACAU GENERALIZADA 


SUBROUTINE CURkEL(MATRIX PETRIA, VECLUR) 

CorsOn GE, N,nES,J 

INTEGER ALFA, GE 

DIFENSLON MATRIAC(NPS,N,02CE),PEIRIACNES, DN USGE), VECTORCN) 
EPSLON=I E-3b | 

JJ=U 

DU 5 JZ=1,h 

LO 5 Ji=1,NPS 

vO 2 J3=0,Gt 

IFUMATRIX(II, 2,93). tU.0)GU TU 2 

IE (JZ Ene 1)GO 10 1 
PeÍnIACUL, JZ, JI)=IIAVECTORCIZ)+EPSLON 

So vu 2 

DO 10 ALFA=O,GF ' 
LFCPETRIXCEATRIXCIL J2,J3) JL ALFA) NE O U)JJ=ZALFA 
CosrInte 

PETRIAX(JL,JZ, JI)D=PLIRIXCMATRIACÇIL  JZ,J3) J2-L, JJ) 
+J3*VECIVUR(JZ)+FPSLON 

CorTiNUE 

PSAXTPETRIACIL,J2,0) 


sa 


13 


Cao 


DO 13 J3=0,6H : 
PLAXTAMAXICEZAS FETRIXCII,J2,J3)) 
CONTANUE 

IELPMAX ES OLU)GO TO 5 

DO 5 J3=0,GÊ 

IFCPMAX EG. PETRIXCJI,)2,93))GU TU 5 
PETRIA(U1,02,93)=0,0 

CONTINUE 

HETURN 

ENU 

CALCULANDO à PALAVRA DECODIE ICELACSEQUENCIA SUBREVIVENTE) 


SuUZKROUTINE SURVIVLFATRIX, PELRIA,hAT) 
CLMMON GE Ma NES, 

INTEGER ALEA,GE, HAT(N) MATRIXCNPS, 1, OSGE) 
DISENSION PETRIXCNPS,Nn, 0IGÊ) 

Inb=1 

DO 2 Kx=N,1,-1 

DO 1 ALFA=O6,6F 
IECPETRIXCIND RX ALFA) NESUCO)HAICAK)=ALEA 
CONTIRUL 

IND=MATRIXCIND, AX BATCKA)) 

CONTINUE 

RETUKN 

LNU 
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APÊNDICE B 


Seja um código de bloco linear binário C(n,kpd) 
com m = 2* palavras códigos, denotadas S; do = 0,l,2pcc els AM 
palavras são transmitidas pela fonte equiprovavelmente atravês de 
um canal perturbado com ruído aditivo branco gaussiano com densida- 


de espectral de potência unilateral No/2. Então 


dicas * Pe leo] =1- P [c leo] (B-1) 


Como & considerado que a ênupla transmitida foi ' 


co! à palavra toda zero, a decodificação por máximo de verossimi - 


lhança serã correta se e somente se 


E cp tn $ps O para cada i=1,2,3,...,m-1 .(B-2) 
t=0 À 


Se os sinais são transmitidos atravês de pulsos '! 


liga-desliga com amplitude +/E Volts ou 0 Volts, então 


p= 
RA (r,-VE/2)< 0 Ve ec ec fe. (B-3) 


Definindo eventos A; para í = 2paseghel; 


E n-1 n-1 
BE Fa Cie n,<v/B/2 E (B-4) 


tem-se que (B-1) pode ser reescrita como (B-5) usando (B-2), 


m=1 
Paioco =1-P [ A, (B-5) 


=L57- 


n-1 
Por outro lado, W(c.) = E c,, êo peso da pala 
ii too É = 
ESA AR E R HA EN, é a relação sinal/ruído do canal. 


Sejam variáveis aleatórias ing, com distri - 


nº Do >n,- Na. (B-6) 


os eventos A, podem ser agora expressos sob a forma 


= i = aa -— — 
A, (N,< YR Yp Tc;)/2) i Td em-1 (B=7) 
m-1 
As variáveis aleatórias (N,) são conjuntamen 


te gaussianas, com matriz de correlação A= (A, 3) onde À, 3"B(N,. N, E; 
portanto a correlação Má serã 
n=-1 a 
PERO dad a 


Ass E (B-8) 
ei vite,).METE,T.Nç/2 


Usando o fato que (np) são variáveis aleatô- 
z 


ê a E 
rias com variância np = N,/21 tem-se que 


E qn2) 
ê c,9º., Ein 
1: aa ddr + Mess (B-9) 
e Yate,T «ATE T «No/2 fale ;) file) 
Notando que | |c; || = E da = (W(c,)) A se- 


gue-se que 


C..l. cao A 
5 =" 1. = cos SpeSs (B-10) 
le; HE He | 
a 
Deve ser observado que 0 < A+3 <l,e 45 
<=> €, ai Ss" ademais As = 1. 


Desta maneira (B-5) pode ser reescrita como sen- 


do 
E — o at? Ex 
Potoco — Em Nasal ( RYpr lc) 7/2, (BY (cq=1) 72) (B-11) 
12 m-1 
Rs é a função distribuição correspondente a 


m-1 variáveis aleatórias conjuntamente Gaussianas' com matriz de co- 


variância A = Ajs)- 


No cálculo da probabilidade de erro por bloco pa 
ra o código (15,11,3), por exemplo, & necessário a avaliação da fun- 


ção distribuição conjunta de 2047 variáveis aleatórias dependentes. 


o cia 


APÊNDICE C 


LEMA: O algoritmo de decodificação por distância mínima generalizada 


escolhe a palavra código que maximiza a confiabilidade média dos sím 


bolos. 


PROVA: A partir das amostras recebidas do canal, dois vetores são de 


terminados, 


1 
E = (ryeser) RO O SPP 1a 
palavra recebida vetor de confiabilidades 
Aqui a é o coeficiente de confiança condicional ! 


de Kiefer associado à decisão r;- 
Se uma palavra código £ é assumida com tendo sido ! 


transmitida, a confiabilidade associada aos seus dígitos & expressa" 


por 
(i) 
Rg se a (rf) =0 
(i) pa 
1-Rs se A (r,1Ê)) dh 


Deste modo, a confiabilidade média dos símbolos da pa- 


lavra, admitindo que a palavra transmitida foi £, é 


Res Late part), + qr-a tr, E Rb 
s nodes) TRE s RESESC Er 


Um algorítmo de decodificação que maximize a confiabi- 


lidade média das decisões deve escolher uma palavra f de acordo com 


n 
Ee = (i) 
Ho ES test I=Rg +41 a (r,1£,))Rg 
n : n - 
(i) (i) 
MARE da RW - x A (x, 1f )12Rç 1) 


fec à=l * i=1 


[60 = 


Lembrando que p tê =— RA - 1, obtém-se 
ES (1) 
Max - E dg (rx 1f,)p 
fec  i=l + 
n (i) 
o que equivale a Min x dy (x 1f,)p 
feC- dis 


A seguir são apresentados dois exemplos ilustrativos. 
A fonte, o canal e código são admitidos os mesmos do exemplo da pági 
na 73. 
Exemplo 1 - Admitindo que as amostras foram(-0.7,-0.2,-0.3,+0.7,+1.1) 
tem-se que 
E =(0,0,0,1,1) e R$ =(.80,.60,.65,.80,.90). 


Dai pode ser encontrado a! dis 


(1040273568) o 

A palavra recebida r não é uma palavra código. Com 
relação a duas palavras código iéS (1,0,0,1,1) e E ESA (E SR Da) Le LD) 
como exemplo, qual delas seria preferível? Aplicando-se o algorítmo, 


a palavra f, seria escolhida: 


> 
Is 
x 
HH 
N 


Com relação aos coeficientes de confiança condicio - 


nal 


) 
) 


sc 
é 

( 
És, R, 


1] 


.63 


(.20,.60,.65,.80,.90) Média RÚD = 3,15/5 


] 
H 


(.80,.40,.35,.80,.90) Média RÉ 3425/85 0H 


Exemplo 2 - Considerando um segundo caso onde as amostras recebidas 


são (-0.7,-0.3,-0.3,+0.7,+1.1), tem-se 


rc = (0,0,0,1,1) e EN e (:80,:65,.65,.80,. 90], 
3056 


portanto p = (bre dge .8) 


=16l= 


não hã preferência. 


= (.20,.65,.65,.80,980) Média R = 3.2/5 .64 


= (.80,.357435,080,.90) Média R5 =3.2/5 = .64 


É realmente indiferente para o valor da confiabilida 
de média dos símbolos se o algorítmo escolhe £, ou £, como sendo a 


palavra decodificada. 


5] 


[6] 


[9] 
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